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^--s ^ This is a review on recent developments of the continuum discretized coupled-channels 

^sj . method (CDCC) and its applications to nuclear physics, cosmology and astrophysics, and 

nuclear engineering. The theoretical foundation of CDCC is shown, and a microscopic reac- 
i2 ' tion theory for nucleus-nucleus scattering is constructed as an underlying theory of CDCC. 

CDCC is then extended to treat Coulomb breakup and four-body breakup. We also propose 

a new theory that makes CDCC applicable to inclusive reactions. 
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§1. Introduction 

Nuclear reaction is one of fundamental reactions in Nature. Understanding of 
nuclear reaction is important not only in nuclear physics but also in cosmology 
and astrophysics. The study of nuclear reaction is necessary to understand the 
mechanism of nucleosynthesis in the early universe, stellar evolution, and supernova 
explosions, and to eventually know the origin of elements in the universe. Ap- 
plication of the study to nuclear engineering is another important direction. The 
Li((i, n) reaction, for example, is regarded as one of the most promising reactions 
to produce intense neutron beams at the international fusion material irradiation 
facility (IFMIF). The study of the Li((i, n) reaction is indispensable for the engineer- 
ing design of accelerator-driven neutron sources. Theoretical calculations play an 
important role there, since experimental data are not available systematically. 

One of the most important current subjects in nuclear physics is to elucidate 
unstable nuclei. This subject is interesting also from the viewpoint of the origin 
of elements through nucleosynthesis. Unstable nuclei are considered to have exotic 
properties such as the halo structure!'''^™ and the loss of magicity in the "island 
of inversion",!^ i.e., the region of unstable nuclei from ^'^Ne to ^"^Mg. In the island 
of inversion, the first-excited states have low excitation energies and large B(E2) 
values.'^''^"'^''^' This indicates that the N = 20 neutron magic number is not valid 
anymore. These novel quantum properties have inspired extensive experimental and 
theoretical studies. 

Important experimental tools of investigating properties of unstable nuclei are 
the interaction cross section ai and the nucleon-removal cross section a-N^^^^^ 
The experimental exploration of halo nuclei is moving from lighter nuclei such as He 
and C isotopes to relatively heavier nuclei such as Ne isotopes. Very lately ai was 
measured by Takechi and collaboratorsli^ for ^*~^^Ne located near or in the island 
of inversion. Furthermore, a halo structure of ^^Ne was reported by Nakamura and 
his group^ with the experiment on (J-n- The ^^Ne nucleus resides in the island of 
inversion and is the heaviest halo nucleus suggested experimentally in the present 
stage. 

The elucidation of unstable nuclei can be accomplished with high-accuracy mea- 
surements of the scattering of unstable nuclei and accurate analyses of the measure- 
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merits with reliable reaction theories. The scattering of unstable nuclei have two 
features. First, the projectile is fragile and thereby the projectile breakup becomes 
important. Second, measurements of the elastic scattering are not easy because of 
weak intensity of the secondary beam, and consequently, there is no reliable phe- 
nomenological optical potential. It is then important to construct a microscopic 
reaction theory that can treat the projectile breakup and does not need phenomeno- 
logical optical potentials. This construction is precisely a goal of the nuclear reaction 
theory. 

A pioneering work on the microscopic description of nucleon-nucleus scatter- 
ing was done by Watson,'i2li which was reformulated by Kerman, McManus, and 
Thaler^^' as series expansion in terms of an underlying nucleon-nucleon (NN) t ma- 
trix. Another important microscopic model particularly to treat inclusive reactions 
is the Glauber model.--' The Glauber model is based on the eikonal and the adia- 
batic approximation. It is well known that the adiabatic approximation makes the 
breakup cross section diverge when the Coulomb breakup is included. For this rea- 
son the Glauber model has been applied mainly for lighter targets in which the role 
of the Coulomb interaction is negligible; see, e.g., Refs . fin]) M ^ HI) . [18]) .fTOl) M^ .[H ]) . 
Recently, a way of making Coulomb corrections to the calculations has been pro- 
posed,'^''^ which works when Coulomb breakup is not so strong. 

The continuum discretized coupled-channels method (CDCC)^''^ is an accu- 
rate method of treating the projectile breakup in exclusive reactions such as elastic 
scattering, elastic-breakup reactions, and transfer reactions. CDCC has succeeded in 
reproducing experimental data on the scattering of both stable and unstable nuclei; 
see, e.g., Refs. |26|) , [27|) , l28l) , l29^, '30^,1^1) ,1321), 1331) ,1^^, 1351^ 

i3]),|MD,ll5]),[l6D,|171),|l8),49),50),51),52),53),54),55),56),57),58),59) and references 
therein. Another reliable method of treating the projectile breakup in exclusive reac- 
tions is the dynamical eikonal approximation (DEA).'S2J'ED The DEA is accurate at 
intermediate and high incident energies where the eikonal approximation is reliable. 

CDCC was first proposed as a method of treating the three-body system com- 
posed of a target and two fragments of a projectile. The original version of CDCC is 
called three-body CDCC. In the review articles IMjl and[25|) on three-body CDCC, 
the method treated nuclear breakup but not Coulomb breakup, and the phenomeno- 
logical optical potentials were mainly used as the interactions between the target and 
the two fragments of the projectile. For the scattering of ^'"^Li and ^^C, there was an 
attempt to obtain the interactions microscopically. However, some phenomenological 
treatments of the imaginary parts of the interactions were necessary in that period. 
Three-body CDCC thus had some limitations. After the review articles, CDCC has 
been developed in several aspects, and applied to various studies in wide research 
fields. The main purpose of the present article is to review these developments and 
applications of CDCC. 

One of the most important developments is the construction of the microscopic 
reaction theory for nucleus-nucleus scatteringS^i based on the multiple scattering 
theory. This gives a foundation of CDCC and other reaction calculations that adopt 
microscopic nucleus-nucleus interactions based on the NN t ox g matrix. The micro- 
scopic reaction theory is reduced to the double-folding model when both nuclear and 
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Coulomb breakup are weak. Then it allows us to practically obtain a microscopic 
optical potential. As mentioned above, CDCC is designed to accurately treat the 
breakup of a projectile into some fragments. The double-folding model can be used 
the input potentials of CDCC, i.e., the microscopic optical potential between the 
target and each constituent of the projectile. 

In ^ we recapitulate the original version of CDCC, i.e., three-body CDCC, and 
show its theoretical foundation, with clarifying the relation between the Faddeev 
method and CDCC. In ^ we present the microscopic reaction theory for nucleus- 
nucleus scattering. We also discuss the validity of the Glauber model and examine the 
localization prescription by Brieva and RooU^ for the nonlocal microscopic optical 
potentials. The double- folding model is then applied to the scattering of stable 
nuclei and unstable Ne isotopes at intermediate energies, in which breakup effects 
are expected to be small. We propose in ^ eikonal-CDCC (E-CDCC) that treats 
Coulomb-dominated breakup processes very accurately and efficiently. Inclusion of 
dynamical relativistic effects in breakup reactions is also accomplished. In ^ we 
extend three-body CDCC to treat breakup processes of a projectile that has a three- 
body structure, i.e., four-body CDCC. We describe also how to smooth the discrete 
breakup cross sections obtained by CDCC. The eikonal reaction theory (ERT) is 
proposed in ^ which makes CDCC applicable to inclusive breakup reactions. CDCC 
and the extended theories are applied to the scattering of unstable nuclei in ^ to 
reactions essential in cosmology and astrophysics in f|8l and to reactions concerning 
nuclear engineering in ^ Finally, we give summary in ^101 

§2. Theoretical foundation of CDCC 

In this section, we recapitulate the original version of CDCC,'^''^ i.e., three- 
body CDCC, and show the theoretical foundation^ ''23 'ESI in order to understand 
what CDCC is. 

2.1. Three-body CDCC 

We consider the projectile that is composed of two particles b and c exactly or 
approximately. In this case, the scattering of the projectile on a target (A) can be 
described by the A-|-b-|-c three-body system. Figured] is an illustration of the three- 
body system and coordinates among the three particles. The three-body scattering 
is governed by the Schrodinger equation 

[H-E]'F = (2-1) 

with the Hamiltonian 

H = Kr + KR + v{r) + U^{r^) + Uc{r,), (2-2) 

where Kr and K^ are the kinetic energy operators associated with r and R, re- 
spectively, and v{r) is the interaction between b and c, while C/b (Uc) is an optical 
potential of the scattering of b (c) on A. 

In CDCC, the total wave function ^Z' is expanded by the complete set of eigen- 
functions of the Hamiltonian h = Kr+v{r) of the b-|-c subsystem. The eigenfunctions 
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Fig. 1. Illustration of the A+b+c three-body system. Vector r^ (7 =b or c) is the relative coordi- 
nate between 7 and A, whereas R and r form a set of Jacobi coordinates. 

are composed of bound and continuum states. The continuum states are character- 
ized by the orbital angular momentum £ and the linear momentum k of the b-|-c 
subsystem, and they are truncated as 

The (.- and /c-truncations are the primary approximation in CDCC. After making 
the truncations, we further discretize the /c-continuum. The model space V' thus 
constructed is described by 

N 

v = Y.\(t^^)m, (24) 

i=0 

where the (j)i represent the bound and discretized-continuum states of h, and A'^ is 
the number of the </>«. The total wave function ^ is hence approximated into 

N 

^ « -^CDCC = V'^ = Y, Mr)Xi{R), (2-5) 

i=0 

where the coefficients Xj(-f^) of the expansion describe motions between the projectile 
in the states (pi and A. The approximate total wave function •Z/cdcc is obtained by 
solving the three-body Schrodinger equation ()2-ip in the model space V' 

V'[H-E]V'^CDCC = 0. (2-6) 

This leads to a set of coupled-channel equations for Xj(ii), 

Af 

[E-Kr- ei]xi{r) = J2i^i\U\4>j)Xj{R), (2-7) 

j 

where U = ?7b(rb) + f^c(rc) and the e, are energies of the b-|-c subsystem in the 
0j. Solving the CDCC equations (|2-7p with the ordinary boundary condition, one 
can get the 5-matrix elements of the elastic scattering and the projectile-breakup 
reactions .'^''^ 
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As for the discretization prescription, three methods have been proposed so 
far: i) the average (Av) method, ES-Ea^ES ii) the midpoint (Mid) method,E3'E3'[6i 
and hi) the pseudo-state (PS) method P 'EaJ -ES, 113 ,IMJ, Hi ,123,153 j^ the Av and Mid 
methods, the A:-continuum is divided into a finite number of bins. In the Av method, 
the continuum states within each bin are averaged into a single state, whereas in the 
Mid method they are represented by a single state at a midpoint of the bin. The 
Av and Mid methods are vahd, since the two methods yield the same converged 
S'-matrix elements as the width A of bins decreases .^3 '129 The Av method is more 
practical than the Mid one, because the former requires less numerical tasks than 
the latter. Since then the Av method has widely been used as a standard way of the 
discretization. 

In the PS method, on the other hand, h is diagonalized in a space spanned by a 
finite number of L^-type basis functions. The resulting eigenstates can be regarded 
as bound states when the eigenenergies are negative and as discretized-continuum 
states when the eigenenergies are positive. The discretized-continuum states are 
thus obtained by a superposition of analytic basis functions in the PS method. This 
makes numerical tasks much easier. We could not derive, however, a smooth breakup 
cross section from the discrete breakup cross section calculated with CDCC. For this 
reason, the PS method was applied to virtual breakup processes in the intermediate 
state of elastic scattering and transfer reactions, but not to breakup reactions them- 
selves. This problem has recently been solved by the smoothing method proposed 
in Refs. [T7|l . HO]) . B8|l . [50|1 . This smoothing method is presented in ^ 

The S'-matrix elements calculated with CDCC depend on the parameters, £maxj 
fcjnax) ^) i-e., the size of the model space P'. This artifact should be removed by 
confirming that the calculated S'-matrix elements converge as the model space is 
expanded. Indeed, the convergence was shown in Refs. [23|) .[25 |) .[67 |) .[68 |) . The next 
question to be addressed is whether the converged S-matrix elements are exact. This 
point is discussed in the next subsection. 

2.2. Relation between CDCC and Faddeev solutions 

As mentioned in ^2.11 CDCC is based on three approximations, the ^-truncation, 
the fc-truncation, and the discretization of the fc-continuum. The ^-truncation is 
most essential among these approximations as shown below. Now we introduce the 
projection operator V that only selects i up to £max- Obviously, V tends to V in the 
limit of large k^^x and small A. The component Vl^ has no asymptotic amplitudes in 
the rearrangement channels. For example, let us consider a simple case of iraax = 0. 
Then, VUV is the average of U over the angle of vector r. After the angle average, 
the potential VUV becomes a function of r and R. Thus, VUV is a three-body 
potential that vanishes at large R and/or large r, so that it does not generate any 
rearrangement channel. 

The insertion of three-body distorting potentials does not change the mathe- 
matical properties of the Faddeev equations. ^'^' Now we consider VUV as such a 
distorting potential in order to obtain the distorted Faddeev equations, 

{E-Kr-KR-v- VUV)^p^ = viiJi, + Vc), (2-8) 
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{E-Kr- 


- Kr - U^)^Ij^ -- 


= (t/b - vu^v)i^A + c/bV'c, 


(2-9) 


{E-Kr- 


-Kr- Uc)'4>c = 


= {U,-VU,V)i>A + U,^lj^,, 


(2-10) 



where V'A, i^h, and ^c satisfy the relation iP" = ipp^ + V'b + V'c- If Eqs. (|2-8p - (|2-10p are 
added, the distorting potential is canceled and the original three-body Schrodinger 
equation ()2-ip is recovered. In an iterative approach to Eqs. ()2-8p - ()2-10p . the zeroth- 
order solution for ipp^ is obtained by setting the right-hand side of Eq. ()2-8p to 
zero. The zeroth-order solution is nothing but !^CDCC- When l^cDCC is inserted in 
Eqs. (|2-9p and ()2-10p . the equations do not generate any disconnected diagram, since 
^P'cDCC has no rearrangement channel in the asymptotic region. Furthermore, the 
subtractions, C/b — VUy^V and Uc — VUcV, sizably weaken couplings of !P'cDCC with 
V'b and V'c- Thus, l^cDCC is a good solution to the three-body Schrodinger equation 
(j2-ip . when ^max is large enough. Very lately, the CDCC solution has directly been 
compared with the Faddeev solution through numerical calculations and it has been 
shown that the two solutions agree very well with each other.l^ 

The discussion mentioned above is made without using the /c-truncation and the 
discretization of A;-continuum. The ^-truncation is thus most essential in CDCC. 

§3. Microscopic reaction theory for nucleus-nucleus scattering 

In this section, we present a microscopic reaction theory for nucleus-nucleus 
scattering, following Ref. [62|) . The validity of the Glauber model is discussed with 
the microscopic reaction theory. When the projectile breakup is weak, the theory 
is reduced to the double-folding model with the effective nucleon-nucleon (NN) in- 
teraction. The microscopic optical potential constructed with the double-folding 
model is applied to measured reaction cross sections for the scattering of stable nu- 
clei and neutron-rich Ne isotopes around 240 MeV/nucleon. Through the analyses, 
it is confirmed that the microscopic optical potential is reliable and hence can be 
used in CDCC as the potentials between a target and fragments of a projectile. This 
microscopic version of CDCC can be applied to many nucleus-nucleus scattering. 

3.1. Microscopic reaction theory 

The most fundamental equation to describe nucleus-nucleus scattering is the 
many-body Schrodinger equation with the realistic NN interaction Vij 

(K + /iP + /iA+ Y^ Vij - E)'F^+^ = 0, (3-1) 

ieP,jeA 

where E is the energy of the total system, K is the kinetic-energy operator for the 
relative motion between a projectile (P) and a target (A), and /ip (/ia) is the internal 
Hamiltonian of P (A). The scattering of P from A can be described with a series of 
multiple scattering in terms of Vij. In the series, one can first take a summation of 
ladder diagrams between the same NN pair. The summation can be described by an 
effective NN interaction tjj in nuclear medium. Taking a resummation of the series 
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in terms of Tjj, one can get the many-body Schrodinger equation with Tip^ 

iK + hp + fiA + Yl ^^i - ^)'^^^^ = ° ' (^"2) 

iePj'eA 

where it has been assumed that the number of pairs (i, j) is much larger than unity. 
We further assumed that the antisymmetrization between incident nucleons in P and 
target nucleons in A can be approximated by using tjj that is properly antisymmet- 
rical with respect to the exchange of the colliding nucleons. The first assumption 
is satisfied well for nucleus-nucleus scattering, and the second one is known to be 
accurate at intermediate and high incident energies .l^'l^ This is an extension of 
the Kerman-McManus-Thaler formalisnnP2' for nucleon-nucleus scattering to nucleus- 
nucleus scattering. 

As mentioned above, r^j describes nucleon-nucleon scattering in nuclear medium. 
A possible simplification of tjj is to replace tjj with the Briickner g'-matrix interac- 
tion, which has been done in many applications; see, e.g., Refs. [63|1 .[73 ]) .17 ^ .175 ]) . [76|1 . 
[77|) .175]) .17^ .150]) . The (7- matrix interaction, however, does not include effects induced 
by finite nucleus, e.g., effects of projectile breakup and target collective excitations, 
because the interaction is evaluated in infinite nuclear matter. In other words, tjj 
is much more complicated than the g(-matrix interaction. Therefore, in general, it is 
not easy to solve the many-body Schrodinger equation p-2p . However, it becomes 
feasible at least in the following cases. 

First, if we consider a) nucleus- nucleus scattering at high incident energies, or, b) 
scattering of lighter projectiles from lighter targets at intermediate incident energies, 
effects induced by finite nucleus such as projectile breakup and target collective 
excitations are small. Then, the double-folding model becomes reliable, with which 
we can analyze the elastic-scattering cross section or the total reaction cross section 
(Tr. We discuss the validity of the double-folding model for the scattering of type b) 
in ^3.41 and ^3.51 It should be noted that in this case the Glauber model also becomes 
reliable, since Coulomb breakup is negligible. The observables to be analyzed with 
the Glauber model are cr or inclusive cross sections such as the one-neutron removal 
cross section. The validity of the Glauber model is discussed in ^3.21 

Second, let us consider a projectile that is a weakly bound system of two nucleus 
b and c. Then the many-body Schrodinger equation ()3-2p is approximated well into 
the three-body Schrodinger equation (j2-ip . in which the potential U^{r.y) (7 =b or 
c) is constructed to reproduce the scattering of 7 on A. When 7 is not a weakly 
bound system, the potential can be constructed microscopically with the double- 
folding model, i.e., by folding the effective NN interaction with the densities of A 
and 7. Thus, we can solve with high accuracy Eq. ()3-2p by three-body CDCC with 
all the input optical potentials obtained microscopically. Note that the microscopic 
potential is nonlocal in general, which makes it very difficult to solve Eq. (|2-ip . It is, 
however, possible to derive a local potential equivalent to the nonlocal potential by 
using the method proposed by Brieva and Rook.l^ The validity of their localization 
method is discussed in 
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3.2. Validity of the Glauber model 

The Glauber model is based on the eikonal approximation for NN scattering 
and the eikonal and adiabatic approximations for nucleus-nucleus scattering. The 
condition for the eikonal approximation to be good for NN collisions in both free 
space and nucleus-nucleus scattering is that 



v(r)/e < 1, A:a>l, 



(3-3) 



where e {k) is the kinetic energy (wave number) of the NN collision, r is the relative 
coordinate between two nucleons and a is the range of the realistic NN interaction 
V. Obviously, this condition is not well satisfied, because v has a strong short-range 
repulsive core; for example, t; ~ 2, 000 MeV at r = for the AV18 force.l^ In fact, 
the eikonal approximation is not good for NN scattering at intermediate energies, as 
shown in the left panel of Fig. [2j To avoid this problem, a slowly-varying function 
such as the Gaussian form has been used as a profile function in the Glauber model.'^ 
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Fig. 2. The on-shell NN scattering amplitude /NN(q) at the laboratory energy -Enn = 150 MeV 
calculated with the bare NN potential AV18 in the left panel and with the JLM g matrix^*' in 
the right panel. The solid (dashed) and dotted (dash-dotted) lines show, respectively, the real 
and imaginary parts of fNN(q) of the exact (eikonal) calculation. 



The use of a slowly-varying profile function and hence of a slowly-varying NN 
interaction can be justified by using the many-body Schrodinger equation (|3-2|) . 
Applying the adiabatic and eikonal approximations to Eq. ()3-2p . one can obtain the 
S matrix of nucleus-nucleus scattering as 



5 = exp 






^^ij '^ij 



(3-4) 



where frei stands for a velocity of P relative to A and Zij is the z-component of 
the relative coordinate Vij between two nucleons. In general, tjj has much milder r 
dependence than the bare NN potential Vij. At high incident energies, for instance. 



n 



is reduced to the t matrix of NN scattering that is a product of Vij and the wave 
operator of the NN scattering. When Vij has a strong repulsive core at small r, the 



u 
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wave operator provides a large suppression there. This leads to the fact that the 
t matrix is a slowly- varying function of rP^ This is also the case with the g matrix. 
The t- or (/-matrix is thus more suitable than Vij as an input of the Glauber model. 
In fact, as shown in the right panel of Fig. [21 the eikonal approximation is quite 
good for NN scattering at intermediate energies, say, 150 MeV, when the g matrix 
proposed by Jeukenne, Lejeune, and Mahaux (JLM/^S is adopted. The use of the 
(/-matrix interaction has another merit in the sense that the effective interaction 
includes, at least in part, nuclear medium effects. 

3.3. Localization of nonlocal microscopic optical potential 

In this subsection, we focus our discussion on nucleon-nucleus scattering for 
simplicity. The multiple scattering theory proposed by Watsonl^ was reformulated 
by Kerman, McManus, and ThalepI^ as series expansion in terms of an underlying 
NN t matrix. The expansion was developed as the spectator expansion,!^ and the 
corresponding first-order optical potential was successful in reproducing experimental 
data in a wide range of incident energies E-m from 65 MeV to 400 MeV.^^' The 
microscopic optical potential thus produced is complex and nonlocal. 

Another way of obtaining the microscopic optical potential is the single-folding 
model with the ^-matrix interaction.E3'E3>E3,ll3,[Zi>E3,IIl-[Z3,[8i xhe ^-matrix in- 
teraction is evaluated in infinite nuclear matter, and the potential between an inci- 
dent nucleon and a target is constructed by folding the g'-matrix interaction with the 
target density by using the local-density approximation. This approach also is highly 
successful in reproducing nucleon-nucleus elastic scattering data^'l^ for 40 MeV 
< -Ein < 800 MeV from light to heavy targets. The optical potential constructed in 
this approach is, again, complex and nonlocal. 
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Fig. 3. The differential cross sections of the proton elastic scattering from ^"Zr at (a) 65 MeV 
and (b) 800 MeV. The lower (upper) horizontal scale shows the transferred wave number q (the 
scattering angle 6cm)- The solid curves represent the results of the exact calculation, whereas 
the dashed lines correspond to those of the equivalent local potential. Experimental data are 
taken from Refs. IM|) and l57)l . 
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Thus, there have been several works on microscopic optical potentials. The use 
of a nonlocal optical potential is, however, not practical in many applications. This 
is true also in CDCC calculations, since the nonlocality makes it difficult to solve 
the CDCC equation ()2-7p . In Ref. [63]), Brieva and Rook proposed an approximate 
form of the equivalent local potential. This Brieva-Rook (BR) localization has then 
been commonly used in many studies ,'^''^''^ the accuracy of which has not yet 
been examined numerically. 

In Ref. |88]) , we tested the validity of the BR localization for the proton elastic 
scattering from ^*^Zr in a wide range of incident energies from 65 to 800 MeV. As 
shown in Fig. [3l the BR localization works very well for (? < 2 fm~ at 65 MeV and 
for Q' < 4 fm~ at 800 MeV. The BR localization is thus accurate in a wide range of 
65 < E^in < 800 MeV, unless q is large. 





Fig. 4. The real part of the exact and approximate wave functions for the p+^"Zr scattering at 
the grazing angular momentum. See the text for details. 



In Fig. m the real part of the exact wave function (solid line) of the Schrodinger 
equation with the nonlocal microscopic potential is compared with that of the ap- 
proximate wave function (dashed line) calculated with the equivalent local potential 
for the p+^*^Zr scattering at the grazing angular momentum. For both E-^^ = 65 MeV 
(left panel) and 185 MeV (right panel), the two results agree very well with each 
other. The exact wave function is thus not suppressed at small R by the nonlocality 
of the optical potential, i.e., the Perrey suppression factor®^' is not necessary in the 
BR-type localization. 

3.4. Application of the double-folding model to the scattering of stable nuclei 

In this subsection, we consider the scattering of stable nuclei from lighter targets 
at intermediate energies. Then, projectile breakup and collective target excitations 
are considered to be small, and consequently, the double-folding model becomes 
reliable. In the double-folding model, the potential U is obtained as a sum of the 
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direct part U^^ and the exchange part U^^^^SiMi 



C/^^(i?) 






Pp('^p)Pa(^a)5°v^(s; Pt,y)drpdrA, 



(3-5) 






X ^{s\ p^u) exp [-zK(i?) • s/M]drpdrx 



(3-6) 



with s = rp — rA + -R- The relative coordinate between P and A is denoted by R, 
and rp (j^a) is the coordinate of a nucleon in P (A) from the center-of-mass (cm.) 
of P (A); /i (z^) stands for the z-component of the isospin of a nucleon in P (A). The 
original form of U is nonlocal, but it has been localized in Eq. ()3-6p with the BR 
localization, where M = A-pAx/{A-p + Aj^) with ^p {Ax) the mass number of P 
(A) . The direct part g^^ and the exchange part g^^ of the g matrix are assumed to 



depend on the local density 



Pij^u = p'^{rp - s/2) + p\{rp^ + s/2) 



(3-7) 



at the midpoint of the interacting nucleon pair. As for the projectile density /9p and 
the target density p\, we use the phenomenological proton-density® deduced from 
electron scattering, and assume that the neutron density has the same geometry 
as the corresponding proton one, since the root-mean-square radii of proton and 
neutron agree with each other with more than 99% accuracy in the Hartree-Fock 
(HF) calculation. 
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Fig. 5. Angular distribution of the ^^C+^^C elastic cross section at 135 MeV/nucleon (left panel) 
and reaction cross sections for ^^C scattering from ^^C, ^"Ne, ^^Na, and ^^Al (right panel). In 
each panel, the solid (dashed) line shows the result with gMB (tLp). The dotted line in the right 
panel represents the result with jmb multiplied by 0.978. The experimental data are taken from 
Ref.lSll) (left panel) and Refs.[ll]),[95l),[96l),[93 (right panel). 

We show in the left panel of Fig. [5]the angular distribution of the ^^C-|-^^C elastic 
scattering at 135 MeV/nucleon. The folding model calculation with the Melbourne 
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g matrix (^mb) reproduces well the data,—' whereas that with the Love-Franey t- 
niatrix interaction (tLp)^ does not. The medium effect is thus important, and the 
double-folding model with ^mb is found to be quite reliable at intermediate incident 
energies. The right panel shows reaction cross sections ctr for ^^C scattering from 
^2C, 20Ne, 23Na, and ^^Al targets at 250.8 MeV/nucleon. The solid and dotted 
lines represent the results of folding model calculations with gyiB before and after 
the normalization oi F = 0.978, respectively. Before the normalization procedure, 
the results (solid line) slightly overestimate the mean values of the experimental 
datgfii''9^''^'5''S3 for A = 20-27. After the normahzation, the results (dotted line) 
agree with the mean values for all the targets. The dashed line corresponds to the 
results of double-folding model calculations with ^lf and no normalization. One sees 
clearly that the medium effect on cjr is also significant at this energy. 

3.5. Application of the double-folding model to reaction cross sections for Ne isotopes 

In this subsection, the double-folding model is applied to the scattering of Ne iso- 
topes from a ^^C target at 240 MeV/nucleon. The projectile densities are constructed 
by either (I) antisymmetrized molecular dynamics (AMD)f^ with the Gogny DIS 
inter actiori22Ji.llDQli or (II) the deformed Woods-Saxon (DWS) model with the defor- 
mation evaluated by AMD. 
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Fig. 6. Reaction cross sections for the scattering of Ne isotopes from a '^^C target at 
240 MeV/nucleon. The experimental data are taken from Refs. Ilip and I96p . The dotted 
line represents the result of the DWS model (method II), whereas the solid line corresponds 
to the result of AMD calculations (method I). The closed square represents the result of the 
AMD-RGM calculation. 

Figure [6] represents the aji for the scattering of Ne isotopes from a ^^C target 
at 240 MeV/nucleon. The result of the double-folding model with the AMD density 
(method I) shown by the solid line reproduces well the datspJi except for ^^Ne. It 
turns out that deformation effects give a significant increase in the ctr and thus very 
important for the reproduction of the data; even no bound state is obtained with a 
spherical HF calculation for j4 > 30. 

The underestimation of method I for '^^Ne comes from the inaccuracy of the AMD 
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density in its tail region. To remove this shortcoming, we describe the ground state 
of '^^Ne by a ^'^Ne+n cluster model including several configurations of the valence 
neutron as well as excited states of ^'^Ne below 10 MeV. We adopt the resonating 
group method (RGM) to include antisymmetrization of nucleons explicitlyP^ The 
result of this AMD-RGM calculation is shown by the square symbol, which repro- 
duces well the experimental data-^^ with no adjustable parameter. It is found that 
the tail correction gives an increase in aji by 35 mb. From these findings, we conclude 
that ^^Ne is a halo nucleus with large deformation. 

Although the AMD-RGM calculation is highly successful, it is quite time con- 
suming. Thus, the DWS model (method II) becomes an important alternative. As 
shown by the dotted line, method II simulates well the results of method I for A ^ 31 
and the result of AMD-RGM for A = 31. An important advantage of method II is 
that the tail correction mentioned above is not necessary. This suggests that the 
DWS model is a handy and promising way of simulating the AMD model with the 
tail correction. The difference between method I and method II for ^^~"^^Ne may 
show the tail correction to method I. 

Further investigations on effects of breakup, dynamical deformation, and reori- 
entation were made in Refs. I102p and I103| ). The findings obtained on a^ for the 
30,31 ]\jg_j_ 12 Q scattering at 240 MeV/nucleon are summarized as follows. 

1) The breakup effects of ^^Ne to ^^Ne+n is 0.7%. 

2) Coupled-channel effects due to the deformation, i.e., couplings to the rotational 
states, of ^*^Ne is less than 1%. 

3) Difference between cjr and the interaction cross section cti, i.e., the total inelas- 
tic cross section, is 0.2% for ^"^Ne. 

4) The reorientation effect due to the intrinsic spin of ^^Ne(3/2^) is 0.2%. 
Thus, these so-called higher-order effects on itr, are found to be very small for 

the present reaction systems. It should be noted, however, that an elastic breakup 
or inelastic excitation to a rotational state itself is an interesting subject of nuclear 
physics. Then, CDCC or ordinary coupled-channel calculation for these processes 
will be necessary. Nevertheless, as mentioned in ^3.H the microscopic reaction theory 
constructed is indispensable, because it gives a theoretical foundation as well as input 
potentials to these reaction calculations. 

§4. Proposal of eikonal-CDCC (E-CDCC) as a method of treating 

Coulomb breakup 

For some decades. Coulomb breakup (or Coulomb excitation) processes have 
been regarded as one of the most effective tools to investigate nuclear proper- 
|-jggfT0l|i xhe standard theory to describe such processes is the virtual photon theory 
(VPT) ,1^2211 'IIQSI or sometimes called equivalent photon method, which is based on the 
picture that a nucleus moves along with its classical trajectory, and dissociates by 
the electromagnetic field constructed by the target nucleus. Although the VPT has 
been highly successful in many cases, various higher-order processes that are not in- 
cluded in the VPT can play some important roles, i.e., nuclear breakup, interference 
between nuclear and Coulomb amplitudes, and multistep breakup processes mainly 
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due to strong continuum-continuum couplings, and so forth. Quantitative evaluation 
of these contributions with more sophisticated quantum mechanical models is very 
important for genuine understanding of the "Coulomb" breakup processes. 

The first application of CDCC to Coulomb breakup was done by Hirabayashi and 
SakuragPS for 208pb(6Li,a+d) at 156 MeV. Later CDCC analyses of 5SNi(SB/Be+p) 
at 26 MeV^*^' and 208pb(SB/Be+p) at 44 and 83 MeV/nucleoiPJ were performed. 
The most essential characteristic of these CDCC calculations is that the number of 
partial waves, Lmax, between the projectile (P) and the target (A) is extremely large; 
typically, Lmax = 15,000 was taken in Ref . [37|1 . This is due to the very long-range 
Coulomb coupling potentials, the dipole components in particular. These calcu- 
lations definitely show the applicability of CDCC to Coulomb breakup processes. 
However, inclusion of such large numbers of partial waves is time-consuming and 
requires very careful numerical treatment of the coupled-channel calculations. 

In Ref. B5|l we proposed a new framework for accurately and efficiently describ- 
ing Coulomb breakup processes, i.e., eikonal-CDCC (E-CDCC). Let us consider a 
reaction system of P and A, in which P consisting of a core nucleus (c) and a valence 
particle (b). E-CDCC uses the Coulomb-eikonal approximation to the scattering 
waves between P and A. Consequently, the wave function of the three-body system 
is given by 

nR,r) = ^^,(r)e-^(— o)<^«^e(6,z)</>C^(6,z), (4-1) 

c 

where ^c(^) is the internal wave function of P with c the channel indices {i, i, m}; 
i > {i = 0) stands for the ith discretized-continuum (ground) state, and £ and 
m are, respectively, the orbital angular momentum between the constituents (c and 
b) of P and its projection on the z-axis taken to be parallel to the incident beam. 
rriQ is the value of m in the incident channel. We here neglect the internal spins 
of c and b for simplicity; see Ref. ,46) for the inclusion of the intrinsic spins of P 
in the E-CDCC formalism, b is the impact parameter (or transverse coordinate) in 
the collision of P and A, which is defined by 6 = \/x'^ + y'^ with R = (x, y, z), the 
relative coordinate of P from A in the Cartesian representation. Note that ^{R, r) 
properly takes into account the dependence of the wave function on the azimuthal 
angle (pji around the z-axis, which corresponds to the coherent choice of the wave 
function in the terminology of the dynamical eikonal approximation (DEA).'S2J''lSlli 

The use of the Coulomb incident wave (f)^ {b, z) instead of the plane wave 
exp(i^c ■ R) in the eikonal approximation is one of the most important features 
of E-CDCC; Kc is the asymptotic wave-number vector of P in channel c from A. In 
actual calculations we use the following approximate asymptotic form 

<^^.(^' ^) ^ ^^^e^(^'=^+^^i-(^'^^-^'^^)), (4-2) 

which is valid for large values of b, with rjc the Sommerfeld parameter for channel c. 
We further assume 

Vii0^^(6,z) ^ iK,{R)4>%Sb.^)-^ (4-3) 
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with 

where fin is the reduced mass of P and A, and Zpe (Z^e) is the charge of P (A). 

As in the usual eikonal approximation, the wave function tpc(b, z) is assumed to 
be a slowly varying function and its second-order derivative Ajitpcib, z) is negligibly 
small. Then we obtain the following coupled-channel equations 

'-^Kf\z)Uf\z) = Y.€'i^) ^2(-) ^{^ e^(^^'-^^)^ (4-5) 



C 



,ib), 



where TZlJ,{z) = {Kc'R - Kc'zy^-' /{K^R - KczY'''^. The reduced couphng potential 
5"^ j (z) is given by 

52(Z) = (^e|f/c + U^\<Pc')r e^("^-™')'^« - ^^^^cC, (4-6) 

where Uq (U\^) is the sum of the nuclear and Coulomb interactions between c (b) and 
A. The explicit form of the coupling potential is given in Refs. B5|l and Il6jl . Note 
that in Eq. ()4-5p b is relegated to a superscript since it is not a dynamical variable. 
Solving Eq. (|4-5p under the boundary condition 

hm v^('')(z) = <5,o, (4-7) 

2— > — OO 

where denotes the incident channel, one obtains the following form of the eikonal 
scattering amplitude^ 

/E = /Ruth^^^ + ^ E f'Llo YLm-mo{K'J, (4-8) 

where jR"*° is the Rutherford amplitude. The partial scattering amplitude f'^^Q is 
defined by 

/l;cO = ;^rtc Y ^^ l"- "' [■ScO,L - Oco\ , (4-9) 



where 



Sco,L= lim4''"^^(z) (4-10) 



with bc;L = {L + l/2)/Kc. The eikonal Coulomb-phase factor HcL is defined by 

nc,L = exp{2ir,,ln{L + l/2)]. (4-11) 

It should be noted that in the derivation of Eq. ()4-8p , we have used 

j^\-^M.t>n^-iK.h^-,Lecos<t>n^^^ = ^J^,((L+ 1/2)0) « 27ri*^y^i^Y^A/(^,0), 

(4-12) 
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with assuming L ^ 1 and small scattering angles 9. Note also that rye In (L + 1/2) 
in the exponent of TicL is the asymptotic form of the Coulomb phase shift (JLirjc)-, 
i.e., 

TdVc) -^ r/eln (L + 1/2) (for L » 7]^). (4-13) 

Furthermore, the formal solution to Eq. ()4-5p will be used as the starting equation 
of the eikonal reaction theory (ERT) described in ^ 

The quantum mechanical (QM) correction in the scattering amplitude can be 
performed if one replaces /^^^o ^°^ small L, say, L < Lq, with the QM partial 
amplitude obtained by conventional QM CDCC, 



f'Q 

JL;cO 



L+i J+io 

E E 

J=\L-i\Lo=\J- 



2Lo + l 



Att 



{iomoLQO\JmQ){£mL mQ—m\Jmo) 



Mo 



'^LLo'5«o)e*^"^+"^«H-)™-"^°, 



(4-14) 



where J is the total angular momentum of the three-body system and S is the 
scattering matrix. This correction is valid if /j^ ^o ~ /l cO with high precision for 
L > Lq ; this is in fact the definition of Lq . We show how this is satisfied below (see 
Fig. |8]). Thus, we use the following hybrid scattering amplitude 



Ruth; 



fcO = /""^"5cO + 



27r 
iK^ 



7 > JL;cO^Lni-mo[I^c) + 
L<Lc 



2tt 



fE 



L>Lc 



cO 



lLm-™o(^c), (4-15) 



which enables one to perform coupled-channel calculations for Coulomb and nuclear 
breakup processes with high computational speed, freely from numerical difficulties 
due to extremely large values of L, and with the same accuracy as of conventional 
QM CDCC. 
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Fig. 7. Angular distribution of the elastic (left panel) and total breakup (right panel) cross sections 
for the *B+^*Ni scattering at 240 MeV. The solid and dashed lines show, respectively, the results 
of the QM and eikonal calculations. The dotted (dash-dotted) line represents the hybrid result 
with Lc = 200 (400). 



In the left and right panels of Fig. [3 respectively, we show the elastic cross section 
(Rutherford ratio) and the total breakup cross section for the '^B-|-^®Ni scattering at 
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240 MeV (30 MeV/nucleon), as a function of the scattering angle 9 in the center-of- 
mass (cm.) frame. Details of the numerical inputs are given in Ref. HSj); we here 
remark only that Lmax = 4, 000 is taken in the calculation. The solid, dashed, and 
dotted lines represent the results of the QM, eikonal, and hybrid calculations, where 
Lq is taken to be 4,000, 0, and 200, respectively. In the right panel the result of 
the hybrid calculation with Lq = 400 is also shown by the dash-dotted line. The 
hybrid calculations with an appropriate value of Lq , namely, 200 (400) for the elastic 
(breakup) cross section, agree very well with the QM results; the difference is only 
less than 1%. Another interesting finding is that the difference between the QM and 
eikonal calculations is quite small up to even 6 = 8°. This encourages one to use 
the eikonal approximation in analysis of reactions around 30 MeV/nucleon. If very 
high accuracy is required as in the study on determination of the astrophysical factor 
Sn, however, inclusion of the QM correction for lower partial waves is necessary (see 
^8.ip . In any case, use of the eikonal amplitude is very effective to save computational 
time and eliminate numerical difficulties. 

It is well known that a theoretical description of breakup reactions at intermedi- 
ate energies requires a relativistic treatment of the reaction dynamics, though only a 
relativistic modification on the kinematics has been usually included. In Refs. I107|) 
and llOSp . we developed a full coupled-channel calculation including a relativistic 
treatment of not only the kinematics but also the dynamics, based on E-CDCC. 
An essential ingredient of the relativistic coupled-channel calculation is the proper 
treatment of nuclear and Coulomb coupling potentials between P and A. We adopt 
the form of the Coulomb dipole and quadrupole interactions shown in Ref. I109|) . 
which is obtained from a relativistic Lienard-Wiechert potential with so-called far- 
field approximation; the validity of this approximation is confirmed in Ref. 1108] ). As 
for the nuclear potential, the conjecture of Feshbach and ZabekP® is adopted. 

It should be noted that E-CDCC can easily incorporate the relativistic Coulomb 
and nuclear coupling potentials, because it adopts the cylindrical coordinate repre- 
sentation. On the other hand, inclusion of these potentials in conventional CDCC, 
which is based on a partial wave decomposition in the spherical coordinate represen- 
tation, is very complicated. Fortunately, however, the dynamical relativistic correc- 
tions are found to be necessary only for large L, where the scattering processes are 
described well by the eikonal approximation. 

Figure [5] shows the real part of the breakup amplitude for ^B+^'^^Pb at 250 
MeV/nucleon as a function of L. As a breakup channel, we choose the s-wave 6th bin 
state, whose breakup amplitude has the largest value. For more details of numerical 
input, see Ref. ll08| ). The solid line represents the result of nonrelativistic QM CDCC, 
which adopts Eq. ()4-14p for all L and has no dynamical relativistic corrections. The 
dotted and dashed lines are the results of E-CDCC, based on Eq. (|4-8p . with and 
without dynamical relativistic corrections, respectively. One sees from the figure 
that at small L, i.e., L < 300, the dashed and dotted lines agree very well with 
each other, and deviate from the solid line. On the other hand, at large L, the solid 
and dashed lines show a very good agreement, as desired, and differ from the dotted 
line. Thus, using the amplitude obtained by nonrelativistic QM CDCC for small 
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Fig. 8. Real part of breakup amplitude for ®B+^''*Pb at 250 MeV/nucleon as a function of the 
orbital angular momentum L between *B and ^"^Pb. The final state of *B is chosen to be 
the s-wave 6th bin state, and the z-component mo of the spin of *B in the incident channel 
chosen as 1. The solid, dashed, and dotted lines show the results of nonrelativistic QM CDCC, 
nonrelativistic E-CDCC, and relativistic E-CDCC, respectively. 



L and that by relativistic E-CDCC for large L allows one to construct an accurate 
coupled-channel framework that includes dynamical relativistic corrections and QM 
effects, i.e., relativistic CDCC. 
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Fig. 9. Double differential breakup cross section (DDBUX) for **B+^°*Pb at 250 MeV/nucleon 
including dynamical relativistic corrections (left panel), and its difference from the calculation 
without relativistic corrections (right panel). 



In the left panel of Fig. [9] we show the double differential breakup cross section 
(DDBUX), (faBu/{d€d9), for the ^B+^ospb reaction at 250 MeV/nucleon calculated 
with relativistic CDCC based on Eq. (j4-15p . e is the relative energy of the two 
fragments, p and ^Be, of ^B after the breakup, and 9 is the scattering angle of the cm. 
of ®B (the p and ^Be system). The right panel shows the difference of the DDBUX in 
the left panel from that calculated with nonrelativistic CDCC, which also is based on 
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Eq. ()4-15p but including no dynamical relativistic corrections. Note that a relativistic 
treatment of the kinematics is adopted in both calculations. One observes a sizable 
increase in the DDBUX due to relativity at forward angles {6 < 0.2°) and around 
e = 1 MeV, where the DDBUX has quite large values as shown in the left panel. 
The increase is indeed large, i.e., several tens of %, which shows the importance 
of including dynamical relativistic corrections in analysis of Coulomb breakup at 
intermediate energies. 

Thus, we have developed E-CDCC for treating both nuclear and Coulomb breakup 
very accurately and efficiently. The scattering amplitude fco to each channel c is ex- 
pressed as a sum of the partial amplitudes /^.^q with L the "classical" orbital angular 
momentum. The QM corrections in fco can be made by replacing the /^.^o with the 
corresponding CDCC results, and the replacement is necessary only for small L. 
E-CDCC calculations with the QM corrections are much faster than CDCC calcu- 
lations, with keeping the same accuracy as of the latter. In addition to that, E- 
CDCC makes it clear the relation between the eikonal picture and the QM picture. 
Furthermore, E-CDCC can incorporate a relativistic Lienard-Wiechert potential as 
a Coulomb interaction and the dynamical relativistic effects generated. E-CDCC 
has thus become a standard method for describing breakup reactions in which the 
Coulomb breakup contribution is essential. 

§5. Proposal of four-body CDCC as a method of treating breakup of 

three-body projectile 

In this section, we present a new version of CDCC to treat four-body breakup 
reactions, in which the projectile breaks up into three constituents, i.e., four-body 
CDCC. In principle, CDCC is a method of treating many-body reaction systems. 
However, except for three-body reactions that involve two-body projectiles, it is 
difficult to solve the CDCC equation. The main problem of the difficulty is how to 
calculate a set of discretized continuum states for a many-body projectile. 

For the discretization procedure, as mentioned in ^ three methods have been 
proposed so far: the average (Av) method, the midpoint (Mid) method, and the 
pseudostate (PS) method. The Av method has been widely used and successful in 
describing three-body breakup reactions. The Mid method is shown to give the same 
results as of the Av method. ^^''^'^' These two methods are, however, difficult to apply 
to many-body reaction systems, where the projectile is to be described as at least a 
three-body system. The main difficulty lies on obtaining the discretized continuum 
states of such projectiles; both the Av method and the Mid method require the exact 
many-body continuum wave functions of projectile. An alternative to the Av and 
Mid methods is the PS method,S3'S3'BSJ,|53,[53,|53 ^^ which the continuum states 

are replaced by pseudostates obtained by diagonalizing the internal Hamiltonian of 
the projectile in a space spanned by L^-type basis functions. One can adopt the 
transformed harmonic oscillator (THO)f^ or the GaussiaEF^'S^ form as the basis 
functions. The validity of the PS method was confirmed for scattering of two-body 
projectiles by the agreement between CDCC solutions calculated with the Av and 
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PS methodsP'HS 

A great advantage of the PS method is that we do not need exact continuum 
wave functions to obtain discretized continuum states, which is an essential point to 
extend CDCC for many-body reaction systems. CDCC with the PS method based 
on Gaussian basis functions was shown to be highly successful in describing the 
elastic scattering of ^He, a typical three-body {a + n + n) projectile, at not only 
high energies but also low energies near the Coulomb barrierp^I''^ see ^7.31 for the 
details. Thus, we have developed a new framework, four-body CDCC, for describing 
a reaction process of a three-body projectile and a target nucleus, explicitly taking 
account of breakup channels of the four-body reaction system. Later, a four-body 
CDCC calculation with the THO basis was also successfully carried out for low 
energy ^He elastic scattering."^ 

There is another difficulty in describing many-body breakup reactions, i.e., how 
to obtain a smooth breakup cross section with respect to the breakup energy e of 
the many-body system. Because of the discretization of the continuum states of 
the projectile, solving CDCC equations provides breakup S'-matrix elements to the 
discretized-continuum states, which give discrete breakup cross sections. On the 
other hand, the breakup cross section in reality is of course continuous as a function 
of e. Thus, one needs a way of smoothing the discretized breakup S'-matrix elements. 

If we could adopt the Av method, the smoothing procedure was straightforward, 
because of the clear relation between the continuum state and the discretized con- 
tinuum one. In fact, there has been an attempt!^ to apply the Av method to ^He 
breakup processes by directly calculating the three-body continuum state with utiliz- 
ing the hyperradial continuum wave functions.— J Though this alternative four-body 
CDCC worked very well for the elastic scattering of ^He, convergence of the breakup 
cross section da/de was not obtained as indicated in Fig. 5 of Ref. [59]) . 

With the PS discretization, we can smooth the breakup cross section by calculat- 
ing the smoothing function {il)^'{e)\<^.y), where ip^'{e) is the exact continuum wave 
function with the incoming boundary condition and $^ is the 7th discretized contin- 
uum state. Details of this smoothing procedure are shown in Refs. IT7|) and I52p for 
three-body reactions and in Ref. I51|) for four-body reactions. This method is much 
more practical than the aforementioned Av method for the three-body projectile, 
because we need only the overlap of ■0'~-' (e) with the spatially damping (compact) 
wave function <P^. To obtain a convergence of da/de with very high accuracy is, 
however, still difficult at present. 

In this situation, very recently, we proposed a new smoothing procedure^ with 
the complex-scaling method.'^^ The advantage of this method is that we do not 
need the exact continuum wave function at all, i.e., only the overlap between two 
compact wave functions is required. Below we describe this new method for the 
scattering of ^He by a target A, as a typical example of four-body reactions. 

The Schrodinger equation for the four-body system of our interest is given by 

[H-EtotWM) = Q^ (5.I) 

where (-I-) indicates the outgoing boundary condition and -Etot is the total energy of 
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the system. The total Hamiltonian H is defined by 

H = KR + hp + UnA + UnA + UaA, (5-2) 

where the internal Hamiltonian /ip of ^He is given by 

hp=Ky + Kr + Vnn + Vna + Vna (5-3) 

within the a + n + n three-body model adopted. The relative coordinate between 
^He and A is denoted by R and the internal coordinates of ®He by a set of Jacobi 
coordinates, ^ = {y,r). The kinetic energy operator associated with R (^) is repre- 
sented by Kr (K^), Vxx' is a nuclear interaction between x and x', and UxA is the 
sum of the nuclear and Coulomb potentials between x and A. 

In CDCC with the PS method, the reaction process is assumed to take place in 
a model space 

P' = j;i<P^)(^,|, (5-4) 

7 

where the (p^ are the eigenstates obtained by diagonalizing hp by L^-type basis 
functions. Details of the calculation for ^^ are shown later in this section. The 
four-body Schrodinger equation is then solved in the model space 

V'[H-Etot]V'\^i,'^cc)=^- (5-5) 

This so-called model space assumption has already been justified as mentioned in 
^ there is no essential difference between three-body CDCC and four-body CDCC 
in this point. In fact, we have confirmed that elastic and breakup cross sections 
calculated with three-body CDCC or four-body CDCC converged with respect to 
expanding the model space.H3'H3'l3iS 

The exact T-matrix element to a breakup state of the a + n + n system is given 
by 

T,ip,k,P) = (4-)(p,fe)x(-)(P)|f/-<r'V^+^ (5-6) 

where U = UnA + UnA + UaA and f^ejic" is the Coulomb interaction between ^He 
and A. P, p, and k are the asymptotic relative momenta conjugate to the coordinate 
R, y, and r, respectively. These three momenta specify the kinematics of the four 
particles to be detected in measurements. The final-state wave functions, iV'e (p, k)) 
and \xe (-P))) with the incoming boundary condition, are defined by 



.(Coul) 



,(-) 



Tr + U;^^"";" - (i^tot - e)\ \x's->{P)) = 0, (5-7) 

[hp-e]\4-\p,k))=0, (5-8) 

where E'tot — e = (^-P)^/(2/^ij) and e = (%)^/(2/Xy) + (M)^/(2/Ur); the fi are the 
reduced masses regarding the relative coordinate given in the subscript. Inserting 
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the approximate complete set Eq. (|5-4p into Eq. ()5-6p . one can get with high accuracy 
the T-matrix element^ 

T,{p,k,P) « Y.(4'Hp,k)\'P^)T^ (5-9) 

with the CDCC T-matrix element 

T, = {<P,xi-\P,)\U - <r'VScc) (5-10) 

to the 7th discrete breakup state ^^ that has the eigenenergy e-y . Equation (|5-9|) has 
been derived by replacing P with P^ in Xs (-P)- The T^ are obtainable by CDCC, 
but it is quite hard to calculate the smoothing function {ip^ {p,k)\<Py) directly as 
mentioned above. Fortunately, however, one can find a way of circumventing this 
difficulty for the calculation of the double differential breakup cross section (DDBUX) 
d'^a/{ded!^p). 

Using Eq. ()5-9p . one can rewrite the DDBUX as 



^ = ^^^LJdp'dk'6{e-e')\TAp',k',P')\' « ^n{e,np) (5-n) 



dedf2_ 
with the generalized response function 



7,7^0 

where G^' is the propagator given by 



n{e,np) = i^^-im[ Yl r,*(cP,|G(-)|^y)ry], (5-12 



G^'^ = lim . (5-13) 

v->+o e — hp — ir] 

The propagator G^' operates only on spatially damping functions (P^. This makes 
the calculation of (^-y|G(~-'|<Py) feasible. 

We now use the complex-scaling method in which the scaling transformation 
operator G{9) and its inverse are defined by 

{r,y\Cmf) = e'''fire^',ye^'), (5-14) 

{f\C-\9)\r,y) = {e-3*^/(re-^^2/e-^)}^ (5-15) 

Using these operators, one can get 

(#^|G(-)|^y) = {'Py\C-\9)G'--^Ci9)\<Py,), (5-16) 

where 

Gi~^ = lim -^ (5-17) 

"^ v-^+o e-h^p-ir] 

with /i^ = G{9)h-pG~^{9). When -vr < < 0, the scaled propagator (^|gJ'^|^') 
is a damping function of ^ and $/ . Note that positive 9 has been taken so far in 
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the complex-scaling method, but in the present situation it must be negative since 
G*-^) has the incoming boundary condition. The scaled propagator Gq is a L^-type 
functionp^ and is then expanded with L^-type basis functions: 



G 



(-) 



E 



(548) 



where (hf 



.J is the ith eigenstate of hp in a model space spanned by L^-type basis func- 
tions, i.e., (<^f|/ip|</'^,) = ef5ii'. Inserting Eq. ([5061) with Eq. (IFTHIl into Eq. ([502]) 
leads to a desired form 



dedOf 



1 (2vr)V|P j^^Tfff 



h^ Pn 



£ — ey 



(5-19) 



with 



7' 



{4^t\Gie)\cPy)Ty, Tf = j;r;(^,|c-i(e)|0f 



(5-20) 



Equation ()5-19p contains no exact three-body continuum states ipe {k,p). This is 
the most important advantage of this new smoothing procedure. 

It should be noted that the use of the response function in the complex-scaling 
method to obtain smooth breakup spectra has already been done in some previous 
studies .l^^'l^J^ The essential feature of the present study is the extension of the 
method for genuine reaction calculations. As mentioned above, the T-matrix element 
T^ in Eq. (j5-20p is the solution of the CDCC calculation for the scattering of a many- 
body projectile by a target nucleus, which thus contains all degrees of freedom of 
reaction dynamics within the model space adopted. 




Fig. 10. Breakup cross section for the d+^*Ni scattering at 80 MeV. The sohd hne is the result of 
the direct smoothing method and the open circles show the result of the new method. 



We tested first the validity of Eq. ()5-19p for a three-body breakup reaction, 
for which the exact breakup T-matrix element T{k,P) = X^„(V' (^)l^n)T'n is 
obtainable by taking the overlap (^^~-'(fc)|^„) directly.'^ This approach is referred 
to as the direct smoothing method below. As an example, we consider the ^^Ni((i, pn) 
reaction at 80 MeV; see Ref. HT]) for the details of the CDCC calculation. Figure [TU] 
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shows the differential breakup cross section da/de, which is obtained by integrating 
the DDBUX of Eq. (|5-19p over the sohd angle i7p. The new method (open circles) 
yields the same result as the direct smoothing method (solid line). The validity of 
the new method is thus confirmed. 

Next, the new method is applied to the ^^C(^He, ann) reaction at 229.8 MeV. 
As for the interactions Vnn and Vna in /ip, we take the GPTIUSi and KKNlsfHSJ po- 
tentials, respectively. These potentials with a Gaussian form reproduce well exper- 
imental data of low-energy nucleon-nucleon and nucleon-a scattering. The particle 
exchange between valence neutrons and neutrons in a is treated approximately with 
the orthogonality condition model.'^^ 



Hi n2 




a 

c=l 




Hi 

o 



ys 



«2 



ra 



(b 

a 

c = S 



Fig. 11. Jacobi coordinates of tiiree rearrangement channels (c = 1-3) of ''He. 



For the diagonalization of hp and /ip, we adopt the Gaussian expansion method 
(GEM).'ii3 In GEM, the a + n + n system is described by a superposition of three 
channels, each of which corresponds to a different set of Jacobi coordinates {y^,rc), 
as shown in Fig. [TTJ For each channel (c), the radial parts of the internal wave 
functions are expanded by a finite number of Gaussian basis functions 



^jxiVc) = y,"e~(^=/S:')V,(%), ^,,{r,) = rfe-(^^/^"»)V,(/^.J, 



(5-21) 



where A (i) is the orbital angular momentum regarding yc (re)- The range parame- 
ters ijj and fj are taken to lie in geometric progression 



Vj = (ymax/yi)^^' 1)/^---, ri = (rmax/r-i)(* ^)/*" 



(5-22) 



The parameters actually depend on c, but we omitted the dependence in Eq. ()5-22p 
for simplicity. In order to confirm the convergence of the breakup cross section with 



Table I. Gaussian range parameters. 



Set 


c 


Jmax 


yi (fm) 


ymax (fm) 


2 111 ax 


n (fm) 


^max (fm) 


I 


3 


10 


0.1 


10.0 


10 


0.5 


10.0 


1, 2 


10 


0.5 


10.0 


10 


0.5 


10.0 


II 


3 


15 


0.1 


20.0 


15 


0.5 


20.0 


1, 2 


15 


0.5 


20.0 


15 


0.5 


20.0 


in 


3 


20 


0.1 


50.0 


20 


0.5 


50.0 


1,2 


20 


0.5 


50.0 


20 


0.5 


50.0 
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respect to the size of the model space, we prepare the three sets of basis functions 
shown in Table HI For the 0+ and 1~ states, maximum internal angular momenta 



in 



and Amax are both set to unity. For the 2+ states, we take 



At, 



for c = 1 and 2, and £max 
Ref.ESI). 



Ar 



2 for c = 3. For other numerical inputs, see 




4 

£ [MeV] 

Fig. 12. Convergence of the breakup cross sections to the 0^ (a), 1~ (b), and 2^ (c) continua. In 
each panel, the dashed line, the solid line, and the open circles correspond to results of sets I, II, 
and III, respectively. The dotted hne in (b) shows the result when Coulomb breakup processes 
are switched off. 

Figure [12] shows the breakup cross sections da/de to the 0^, 1~, and 2+ continua. 
For all the cross sections, sets II and III yield the same result, whereas the result 
of set I is somewhat different from it. The convergence of CDCC solution with 
respect to expanding the model space is thus obtained with set II. Figure [13] shows 
the breakup cross sections for four values of the complex-scaling angle 9. All the 
breakup cross sections converge at 9 = —14°, when 9 decreases from —6° to —18°. 



(a) 

e=-6deg 
e=-10deg 
e=-14deg 
e=-18deg 




4 
£ [MeV] 




0.8 1.2 

£ [MeV] 



Fig. 13. Dependence of the breakup cross section to the 0^ (a), 1 (b), and 2"*" (c) continua on 
complex-scaling angle 0. 

The present calculation includes nuclear and Coulomb breakup processes. The 
contribution of the latter is indeed large for the breakup to the 1^ continuum, as 
indicated by the dotted line in Fig. [T^Tb). It increases the breakup cross section 
by a factor of 2 from the result without Coulomb breakup. On the other hand, 
the Coulomb breakup effect turns out to be negligible for the 0"*" and 2+ breakup 
spectra, and the latter dominates the breakup cross section. Consequently, the 
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Coulomb breakup effect in ttie present reaction system is not significant for botfi tfie 
breakup reaction and the elastic scattering. For heavier targets, however, Coulomb 
breakup processes become dominant. Four-body CDCC can treat both nuclear and 
Coulomb breakup processes and hence their interference. In m.2>\ we apply four- 
body CDCC to ®He scattering from both light and heavy targets in a wide range of 
incident energies. 

§6. Proposal of the eikonal reaction theory (ERT) that makes CDCC 

applicable to inclusive reactions 

CDCC is highly reliable for describing exclusive reactions but not applicable 
to inclusive cross sections such as a neutron removal cross section. In this section, 
we present an accurate method of treating inclusive reactions within the framework 
of CDCC. According to the method, the nuclear and Coulomb breakup processes 
are consistently treated by CDCC without making the adiabatic approximation to 
the Coulomb interaction, so that the removal cross section calculated never diverges. 
This method is referred to as the eikonal reaction theory (ERT). This section consists 
of a brief review of Refs. [TT9|) and ll20p . and a new application of ERT to two- neutron 
removal reactions. 

6.1. Formulation of ERT 

Let us consider as a projectile (P) the two-body system composed of a core 
nucleus (c) and a valence neutron (n), and hence the three-body (c+n+A) system 
for the scattering of P from a target A. We start with the 5- matrix operator as a 
formal solution to the coupled-channel (E-CDCC) equations ()4-5p in fJU 



exp 



iC / dzO^UO 



(6-1) 



with the interaction 



U = t/^N"=i) + C/^N^^i) + C/iC°"i) , (6-2) 



where Un "'^ is the nuclear part of the optical potential between n and A, and 
Uc ^^ and Uc °^ are, respectively, the nuclear and Coulomb parts of the optical 
potential between c and A. C is the path ordering operator which describes the 
multistep scattering processes accurately. The operator O is defined by 

6 = ^e^^■^ (6-3) 

\Jhv 



where the wave-number operator K is given hy K = y^2^{i{E — h)/h with the in- 
ternal Hamiltonian h of P, and v = hK /fipi. is the velocity operator. In the definition 
of O we assume that there is no Coulomb monopole interaction between P and A 
for simplicity. In the actual calculation, the exponent in Eq. ()4-2p . with Kc and 
T]c replaced with the operator form, is used as an exponent in O; the form of v is 
changed accordingly. 
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In the Glauber model, the adiabatic approximation is made, in which h is re- 
placed with the ground-state energy eo, and hence O^UO and C in Eq. (j6-ip are 
reduced to U/{hvo) and 1, respectively, where vq is the velocity of P in the ground 
state relative to A. At intermediate energies, this treatment is known to be valid 
for short-range nuclear interactions, but not for the long-range Coulomb interac- 
tions. Therefore, we make the adiabatic approximation in the evaluation of only 
O^Un ^'^ O, i.e., we use 



Qt[/^Nucl)Q ^ UJi^'^^^^/ihvo). 



(6-4) 



For the scattering of '^^Ne from a ^ospj^ target at 240 MeV/nucleon, it is shown 
by CDCC that the effect due to the replacement (|6-4|) is small; it is 0.2% for the 
reaction cross section a^i, 1.9% for the elastic-breakup cross section (TeBj 4.1% for 
the one-neutron stripping cross section o"„:STR- Using this replacement, S can be 
separated into the neutron part Sn and the core part 5c, 



with 



Sn = exp 
^c = exp 



O — Jn^c 



^ / dz^i^ucl) 



hvo 

iC 



(izOt(C/iN"=') + t/i^°"'))0 



(6-5) 

(6-6) 

(6-7) 



This is the most important result of ERT. We can derive several kinds of cross 
sections with the product form Eq. ()6-5p . following the formulation of the cross 
sections in the Glauber model™''^ see, e.g., Ref. I120p for the explicit form of 
them. It should be noted that one cannot evaluate 5c directly with Eq. ()6-7p . since 
it includes the operators O and C. However, one may find that 5c is the formal 
solution to the Schrodinger equation 



J^ 
2/iij 



Vl + h + C/(N"=') + [/(C-i) - E 



l^c = 0. 



(6-8) 



when the eikonal approximation is made. We can thus obtain 5c by solving Eq. (j6-8 
with E-CDCC.'^''^ As mentioned above, 5„ is obtained by Eq 



6.2. One-neutron removal reaction of^^Ne 

We apply ERT to the one-neutron removal reactions for the ^^Ne+^^C scat- 
tering at 230 MeV/nucleon and the ^^Ne+^'^^Pb scattering at 234 MeV/nucleon 
with a ^"^Ne+n+A three-body model. The optical potentials for the n-target and 
^'^Ne-target subsystems are obtained by folding the effective nucleon-nucleon inter- 
actioiP^ with one-body nuclear densities. The densities of P and A are constructed 
by the same method as in Ref. I122p . We assume the ground state of ^^Ne to be ei- 
ther the ^°Ne(0+) lp3/2 or the ^°Ne(0+) (g) 0f7/2, with the one-neutron separation 
energy of 0.33 MeV. As for the breakup states, we include s-, p-, d-, f-, and g-waves 
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up to the relative momentum between ^'^Ne and n of 0.8 fm^^. For more detailed 
numerical inputs, see Ref. I119p . 

Table [III presents the integrated elastic-breakup cross section cr^g, the one- 
neutron stripping cross section <T„.gxR) the one- neutron removal cross section 0"_„, 
and the spectroscopic factor S = a^^^/a^\. S calculated with the lp3/2 ground- 
state neutron configuration little depends on the target and less than unity, but that 
with the 0f7/2 configuration does not satisfy these conditions. Therefore, we can 
infer that the major component of the ground state of ^^Ne is ^'^Ne(0"'") ® lp3/2 with 
S ~ 0.69. This value is consistent with the result shown in Sec. 13.51 We adopt this 
configuration in the following. 

Table II. Integrated cross sections and the spectroscopic factor for the ^^Ne-^'^C scattering at 
230 MeV/nucleon and the "^^Ne-^^^Pb scattering at 234 MeV/nucleon. The cross sections are 
presented in unit of mb and the data are taken from Ref. I12|l . 





^^C target 


2DB 


Pb target 




P3/2 h/2 Exp. 


P3/2 


h/2 Exp. 


0"n:STR 


23.3 3.3 
90 29 


799.5 
244 


73.0 (540) 
53 


s 


114 32 79 
0.693 2.47 


1044 
0.682 


126 712 
5.65 



Since the potential between ^'^Ne and n is not well known, we change each of 
the potential parameters by 30% and see how this ambiguity of the potential affects 
the resulting value of S. We obtain S = 0.693 ± 0.133 ± 0.061 for a ^^C target and 
S = 0.682 lb 0.133 lb 0.062 for a ^ospj-, target; the second and third numbers following 
the mean value stand for the theoretical and experimental uncertainties, respectively. 
Thus, S includes a sizable theoretical uncertainty. This situation completely changes 
if we look at the asymptotic normalization coefficient (7,^^2211 {^^ (j = 0.320 ib 0.010 ib 
0.028 fm-i/2 for a ^^C target and C = 0.318 ib 0.008 =b 0.029 fm'^s for a ^ospb 
target. Thus, C has a much smaller theoretical uncertainty than S. This means 
that the one-nucleon removal reaction is quite peripheral. 

6.3. Two-neutron removal reaction of ^He 

ERT is applied to two- neutron removal reactions of ^He from ^^C and 208p|-, 
targets at 240 MeV/nucleon. In this case, the projectile is treated as a three-body 
[a + n + n) system and hence four-body CDCC is used. The optical potentials for the 
n-target and a-target subsystems are calculated by folding the Melbourne nucleon- 
nucleon (/-matrix interaction^ with the densities obtained by the spherical Hartree- 
Fock (HF) calculation with the Gogny DIS interaction.'^ 'l^^ The model space of 
the following calculation is the same as in Ref. I53|) , with which good convergence is 
achieved. 

Table IIIII shows the one- and two-neutron stripping cross sections, cr^. g^R ^-nd 
CTgnSTRi respectively, and the two-neutron removal cross section o"_2„. The present 
framework reproduces well the experimental data^^ with no adjustable parameters. 
Thus, we can clearly see the reliability of ERT for two-neutron removal reactions on 
both light and heavy targets. 
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Table III. Integrated cross sections for two-neutron removal reaction of ^He on ^^C and ^"^Pb 
targets at 240 MeV/nucleon. The cross sections are presented in unit of mb and the experimental 
data are taken from Ref. 11241 ). 





^^C target 


""«Pb target 




Calc. 


Exp. 


Calc. 


Exp. 


""niSTR 
''■2n:STR 


153.4 
29.0 


127 ± 14 
33 ± 23 


353.6 
148.9 


320 ± 90 
180 ± 100 


0--2n 


198.5 


190 ± 18 


1016.6 


1150 ± 90 



6.4. Comparison of ERT and the Glauber model 

In this subsection, the accuracy of the Glauber model is tested for deuteron 
induced reactions at 200 MeV/nucleon. We estimate the relative difference 



5x = [X(CDCC) - X(GL)]/X(CDCC), 



(6-9) 



where X(CDCC) and X(GL) correspond to integrated cross sections calculated by 
CDCC with ERT and the Glauber model, respectively. For the Glauber-model cal- 
culation in which the eikonal and adiabatic approximations are made, the Coulomb 
interaction is neglected. The numerical inputs are given in Ref. 1120) . 
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Fig. 14. Accuracy of the Glauber model for various integrated cross sections. 

Figure [HI shows bx as a function of target mass number A for cteb, <7r, (T^iStr, 
(T_n, the proton stripping ((Tp:STR) and removal ((T_p) cross sections, the total fusion 
cross section (Ttf, and the complete fusion cross section ucf- The Glauber model is 
good for light targets, as expected. For heavier targets where the Coulomb breakup 
is essential, however, the Glauber model is not good for the elastic-breakup and 
proton stripping cross sections. 

§7. Applications of CDCC to scattering of unstable nuclei 



In this section, we review some recent applications of CDCC to scattering of 
unstable nuclei. See the references cited in the following subsections for the details 
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of the formalism, numerical inputs, other results, and further discussion. 

7.1. One-neutron removal reactions of ^®C and ^^C on proton target 

One-neutron removal reactions have played a key role in investigating exotic 
properties of neutron-rich nuclei. For light targets such as Be and C, the so-called 
stripping modeP^''^^ based on the eikonal approximation has been successful in 
describing the removal process. On the other hand, for a proton target, the removal 
process can be interpreted as an elastic breakup, since there is no excited state in 
the target, and also a neutron transfer process producing deuteron hardly occurs. 
Therefore, we should treat accurately the elastic breakup process for the one-neutron 
removal reactions on a proton target, to which the stripping model is not applicable. 
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Fig. 15. Level scheme and transitions of ^*C (left panel) and transverse-momentum distribution of 
^**C in coincidence with the 2.4 MeV 7-ray (right panel). The solid and dashed lines represent 
the CDCC calculations by assuming neutron removal from the £ = and 1 = 2 orbits in ^^C, 
respectively. 

In Ref. I127p . one neutron removal cross sections of ^^C and ^^C on a proton 
target were measured, and CDCC was applied to analyze the data to determine the 
neutron configuration of these nuclei. Figure [15] (right panel) shows the transverse- 
momentum distribution of ^^C after the breakup of ^^C in coincidence with the 
2.4 MeV 7-ray, which corresponds to the population of the 4.0 MeV state of ^^C; see 
the level scheme in the left panel of Fig. [T^l The experimental data are consistent 
with the CDCC calculation for neutron breakup from the Od orbit in ^^C(l/2^). 
From this analysis, we have successfully confirmed the spin-parity I'^ = (2,3)"^ as- 
signment for the 4.0 MeV state in ^^C; see Ref. 1127] ) for more detailed discussion. 

For both reactions induced by ^^C and ^^C, the theoretical cross sections cal- 
culated by CDCC with the shell-model spectroscopic factors obtained by the WBP 



interactiorM^^ explain the most of the relative values, as shown in the Tables of 
Ref. ll27p . As for the absolute values, the measured and calculated results agree very 
well with each other for the breakup of ^^C that has a small neutron separation en- 
ergy Sn- On the other hand, the measured cross section is significantly smaller than 
the calculated one for the breakup of ^^C that has a quite large S„. This is consistent 
with the finding of a systematic analysis of nucleon removal cross sections ,'I2J' though 
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the mechanism of this "quenching" of the shell-model spectroscopic factors is still 
under discussion. Thus, CDCC is shown to be a powerful method for determining 
configurations of valence neutrons in exotic nuclei. 

7.2. Study on the low-lying states in the unbound nucleus ^^Be 

Shell evolution of neutron-rich nuclei is one of the hottest topics of nuclear 
physics. In Ref. 1129] ). we investigated the shell structure of ^^Be, which is unbound 
and has a neutron number A^ = 9. We did the same analysis as in ^7.11 of the 
one-neutron removal cross section of ^^Be by a proton target. 
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Fig. 16. Transverse momentum distributions of the ^^Be-|-n system for (a) 0.25 MeV < E-^d < 
0.75 MeV and (b) 2.0 MeV < E^^i < 2.5 MeV. The s-, p-, and d-wave components are shown 
by the dotted, dashed, and dash-dotted hues, respectively. The solid hues are the sums of the 
three components. 



We show in Fig.[Tn]the transverse momentum distributions of the ^^Be-|-n system 
for 0.25 MeV< ^^ei <0.75 MeV (left panel) and 2.0 MeV< ^^ei <2.5 MeV (right 
panel), where E'rei is the relative energy of ^^Be and n. The dotted, dashed, and dash- 
dotted lines correspond to the s-, p-, and d-wave contributions calculated by CDCC, 
respectively. We adopted a ^^Be-|-n two-body model with each of these relative 
angular momenta as an initial state structure of ^^Be, and then added all the results 
incoherently. The strengths of the individual components as well as the resonance 
parameters assumed were treated as free adjustable parameters to reproduce the 
experimental data; see Ref. I129|) for more details. From this analysis, we confirmed 
the low- lying p-wave resonance at E^e\ = 0.51(1) MeV as well as the d-wave resonance 
at E'rei = 2.39(5) MeV on the broad s-wave distribution with the s-wave scattering 
length Os = —3.4(6) fm in the invariant mass spectrum of ^'^Be. The spin-parity of 
the p-wave resonance state was assigned to I^ = 1/2". This intruder l/2~ state 
indicates the disappearance of the N = 8 shell closure and the shell evolution in the 
vicinity of the neutron drip line. 

7.3. Elastic and breakup cross sections of^He 

Two-neutron halo nuclei such as ^He and ^^Li have exotic properties, i.e., soft 
dipole excitation and a di-neutron correlation. These properties can be investigated 
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via breakup reactions, where the projectile breaks up into three (core+n+n) frag- 
ments. For this purpose, accurate description of the four-body systems is essential, 
and four-body CDCC proposed in ^is one of the most reliable methods for treating 
such four-body breakup processes. 
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Fig. 17. Angular distributions of the elastic differential cross section for the ^He+^^C scattering at 
18 MeV (left panel) and 229.8 MeV (right panel). The sohd (dashed) line shows the result with 
(without) breakup effects. The experimental data in the left and right panels are taken from 
Refs. nnO)l and ll3l |). respectively. 

First, we apphed four-body CDCC to the elastic scattering of 6He.S9S Y'lg- 
ure [T7] shows the results of the four-body CDCC calculation compared with the 
experimental data for the ^He+^^C scattering at 18 MeV^^ and 229.8 MeV.HSJ in 
this analysis, coupling potentials in four-body CDCC are calculated by a double- 
folding model with the density-dependent M3Y interaction-'-^^' that has only the real 
part. The imaginary part of each coupling potential is assumed to have the same 
form as of the real part with a normalization constant Ni to be optimized to repro- 
duce the experimental data. We took A^i = 0.3 and 0.5 for the scattering at 18 MeV 
and 229.8 MeV, respectively. One sees clearly from Fig. [T7] the importance of the 
breakup effects; note that A^i does not change the oscillation pattern. In this work, 
we omitted couplings to 1~ continuum states, because the Coulomb breakup effect 
is negligible in this reaction system. 

In Fig. [TSl angular distributions of the elastic differential cross section for the 
6jjg_^209g^ scattering at 19 MeV and 22.5 MeVf^S'lSS are shown. Since Coulomb 
breakup effects are dominant in this system, we take into account couplings to 1^ 
continuum states. In this calculation, we adopt phenomenological optical potentials 
for the interactions between the ^°^Bi target and the constituents of ^He. One 
sees that the four-body CDCC calculation reproduces well the experimental data. 
Again, the breakup effects are significant. It was found in Ref. [50]) also that three- 
body CDCC, in which a di- neutron plus a structure was assumed for ^He, could not 
reproduce the data. This shows importance of the accurate description of ^He by 
means of the three-body model. 

Next, four-body CDCC with the new smoothing method described in ^is ap- 
plied to analyses of breakup reactions of ^He. In Fig. [T^l the calculated breakup 
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Fig. 18. Same as in Fig.[n]but for the *'He+^°*'Bi scattering at 19 MeV (left panel) and 22.5 MeV 
(right panel). The experimental data are taken from Refs. IT33|) and 1134 () . 
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Fig. 19. Comparison of the breakup cross section calculated by CDCC (solid line) with experimental 
data for (a) the "^He+^^C scattering at 240 MeV/nucleon and (b) the ^'He+^^^Pb scattering at 
240 MeV/nucleon. The dash-dotted, dotted, and dashed lines correspond to the contributions of 
the 0"*", 1~, and 2^ breakup, respectively, and the solid line is the sum of them. The experimental 
data are taken from Ref. I124ll. 



cross section da/de is compared with the experimental data for the ^He+^^C and 
^He+^^'^Pb reactions at 240 MeV/nucleonP^ We found that nuclear breakup was 
dominant in the former, whereas Coulomb breakup to the 1~ continuum was dom- 
inant in the latter. For the ^^C target, the present theoretical result is consistent 
with the experimental data except for the peak of the 2"*" resonance around e = 1 
MeV. For the ^osp^-, target, the present method undershoots the experimental data 
for e > 2 MeV. A possible origin of this undershooting is the contribution of the in- 
elastic breakup reactions not included in the present calculation. It was reported in 
Ref. I135|) that the inelastic breakup effect was not negligible, and the elastic breakup 
cross section calculated with four-body DWBA also underestimated the data. 

One of the most important features of these analyses with CDCC is the treatment 
of the nuclear and Coulomb breakup amplitudes as well as their interference on the 
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same footing. Through a systematic analysis of *B breakup reactions, ^^' we showed 
the importance of the nuclear-Coulomb interference which has been neglected in 
many studies on breakup of unstable nuclei. 

§8. Applications of CDCC to reactions essential in cosmology and 

astrophysics 

In this section, we review our recent applications of CDCC to some reactions 
important for cosmology and astrophysics. Readers are invited to the references 
cited in the following subsections for more detailed description of the calculation 
and complete discussion. 

8.1. Determination of the astrophysical factor Sn for the "^ Be(p,j)^B reaction 

Intensive measurements of *B breakup at intermediate energies^^''^'^^^ have 
been done to indirectly determine the astrophysical factor 5i7(0), the value in the 
zero-energy limit, for the Be(p, 7)^B reaction, which is important for neutrino- 
physics.'I^SJ The results of the indirect measurements of 5*17(0) are, however, signif- 
icantly smaller than the result of the precise direct measurement of ^Be(p, 7)*B.'i^ 
This may cast doubt on the reliability of the indirect measurements of astrophysical 
factors through breakup reactions of unstable nuclei, which are planned to be done 
at forthcoming RI beam facilities such as FAIR at GSI and FRIB at MSU, and at the 
brand-new facility RIBF at RIKEN. Thus, to clarify the reason for this discrepancy 
is a very important subject of nuclear astrophysics. For this purpose we accurately 
analyzed the ^B breakup by ^"^^Pb at 52 MeV/nucleon with full coupled-channel 
calculation, i.e., CDCC, taking account of both nuclear and Coulomb breakup with 
all higher-order processes.!^ 
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Fig. 20. a) Result of the CDCC analysis of the *B breakup cross section at 52 MeV/nucleon.^^ 
b) Summary of the present study on Si7(0) compared with the previous results. See the text 
for details. 

In the left panel of Fig. [501 we show the result of the CDCC analysis of the ^B 
breakup cross section by ^"^^Pb at 52 MeV/nucleon.'^^ The cross section, which is 
integrated over the relative energy en of ^Be and p between 500 and 750 keV, is 
shown as a function of the scattering angle 9 of the center-of-mass of ^B (the "^Be+p 
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system). Numerical inputs for the CDCC calculation are described in Ref. 46) in 
detail. The dashed line is the result of the calculation using a normalized ^B internal 
wave function, whereas the solid line is the result of the x^ fitting to the experimental 
data. The theoretical results are smeared out by using a filtering tabl c^ ^ ^ to take 
account of the experimental resolution and efficiency. Since the filtering table was 
made with assuming the s-state breakup of ^B, which turns out to be valid for 
9 < 4°, we use the data in that region in the x^'fittiiig procedure. Note that 
each horizontal bar put on the data points below 4° does not represent a statistical 
error but it shows the range of 9 in which the breakup cross sections contribute to 
each data point .'^^SI xhe purpose of the x^ fitting is to determine the tail of the 
overlap function /(r) between the ground states of ^Be and ^B, i.e., the asymptotic 
normalization coefficient (ANC) C; the resulting value of C is 0.74 fm~^". 

It is shown in Ref. I123p that £'17(0) for the '^Be{p,j)^B reaction is accurately 
determined if C is obtained from an alternative reaction in which only the tail 
region of the overlap function /(r) has a significant contribution. We use this so- 
called ANC method and 5'i7(0) = 20.9^1^ (theor) ±1.8 (expt) eV b is obtained; the 
theoretical uncertainties are carefully evaluated by changing the numerical inputs as 
described in Ref. [36|) . T his value of 5i7(0) is significantly larger than the previous 
one, 19.0 it 1.8 eV b,'^3SJ which was obtained from the same experimental data as 
used in the present analysis with the VPT. It is also found that if we assume simple 
one-step El transition in our analysis of the *B breakup process as in the VPT 
analysis, 19.0_^'7 eV b is obtained,'^ which agrees very well with the previous 
result. This clearly shows the importance of the accurate description of the breakup 
process, taking account of both nuclear and Coulomb breakup with all higher-order 
processes. 

In the right panel of Fig. [20] we summarize the results of the present analysis of 
the ^B breakup reaction, compared with the results of the precise direct ^Be(p, 7)^B 
measuremen4i32l and the previous indirect ®B breakup measurement.'^^ One sees 
that the agreement between the values of 517(0) obtained from the direct and indi- 
rect measurements is significantly improved. It should be noted that very recently, 
5*17(0) = 20.9 ± 0.7 (theor) ±0.6 (expt) eV b was obtained,!^^ as a compilation 
result of direct measurements, with a small correction to the result of Ref. I139p . 

Thus, we conclude that the indirect ^B breakup measurement indeed enables 
accurate determination of the astrophysical factor 5*17(0) if the breakup reaction is 
analyzed accurately. This conclusion will be important for a future project to extract 
nuclearastrophysics information from breakup reactions of unstable nuclei. 

8.2. Three-body model analysis of subbarrier a transfer reaction 

Recently, indirect measurements of ^^C(a, n)^®0 using subbarrier a transfer re- 
action ^'^C(®Li,d)^ 0(6.356 MeV, 1/2"'') has been performed by Florida State Uni- 
versity groupp^ DWBA was used for the analysis of the transfer cross section. ^Li 
is known to have aa+d two-body structure with a small binding energy of 1.47 MeV. 
Furthermore, the final state of ^^0(6.356 MeV, 1/2"*"), which is denoted by "'^'^O* be- 
low for simplicity, locates just 3 keV below the a-^^C threshold. Therefore, roles of 
the breakup states of ^Li and ^^O* in the a transfer process ^^C(^Li, (i)^^O* should 
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be clarified, to obtain a more proper value of the cross section of ^^C(Q,n)^^0 at 
very low energies. 

In this studyp^J ^^q f^j-g^ calculate with CDCC the elastic cross sections of ^Li- 
^■^C at 3.6 MeV and d-^^0* at 1.1 MeV corresponding to the initial and final channels, 
respectively, of the ^^C(^Li, (i)^^O* reaction. It is found numerically that the breakup 
effects of ^Li are important but those of ^^O* are small. Therefore, we evaluate the 
^'^C(^Li,d)^ O* cross section at 3.6 MeV including the breakup states of ^Li with 
CDCC; the transition of the transfer process is described by Born approximation, 
i.e., we adopt the so-called CDCC-BA framework. In the left panel of Fig. [211 we 
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Fig. 21. (left panel) Cross section of the transfer reaction "C(^Li, d)^''0* at 3.6 MeV. The sohd 
line is the result of the calculation with a normalized a-^'^C wave function and the dashed line is 
the result of the ^ fitting to the experimental data.l^^I} (right panel) Dependence of the cross 
section on the model space of CDCC; see the text for details. 



show the transfer cross section calculated by CDCC-BA with a normalized a-^^C 
wave function (solid line) and that after the x^ fitting (dashed line) to the experi- 
mental data.'I^ Through this analysis, we obtained 1.03 it 0.29 fm~^ for the square 
of the Coulomb-modified ANCP^ of ^ O* with the a-^^C configuration, where the 
second number shows the theoretical and experimental uncertainties together. This 
result is consistent with the previous result 0.89 ±0.23 fm~^ obtained by the DWBA 
analysis.'I^ It is found that the transfer processes through the breakup states of 
^Li are negligible and only the so-called back-couplings from the breakup states to 
the ground state of ^Li are important. This means only a proper description of the 
scattering wave between ^Li and ^'^C in the elastic channel is necessary. Therefore, 
use of a reliable optical potential between ^Li and ^^C, as in the previous DWBA 
analysis, may give a proper value of the ANC including the back-coupling effects 
implicitly. However, this will not always be the case, since usually the optical poten- 
tial is determined only phenomenologically. Furthermore, breakup transfer processes 
might be important in other subbarrier a transfer reactions. The present three-body 
approach, therefore, should be applied systematically to these reactions. 

Before ending this subsection, we discuss the convergence of the subbarrier trans- 
fer cross section with CDCC. In the right panel of Fig. [21] we show the dependence of 
the cross section on the maximum value /cmax of the relative wave number of the d-a 
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continuum adopted in CDCC. The values of fcmax are shown in unit of fm~^ and the 
corresponding values of the d-a relative energy emax are given in the parentheses in 
unit of MeV. One sees clearly that the convergence is very slow and obtained even- 
tually at /cmax = 2.0 fm~^ (cmax = 62.4 MeV). In the usual CDCC calculation, one 
takes only the open channels, in which the relative energy Sj-ei between the projectile 
and the target is positive. The result thus obtained (thin solid line) is, however, siz- 
ably different from the converged one (thick dotted line). Therefore, inclusion of the 
closed channels, in which E'l-ei is negative, is necessary for the accurate description 
of the i3c(6Li,(i)i^O* reaction at 3.6 MeV with CDCC. 

8.3. Nonresonant triple a process at low temperatures 

Understanding of the formation of ^^C, an essential element of life, is one of 
the most important subjects in physics. It is well known that the so-called triple-a 
process, i.e., the sequence of the following two reactions 
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is the path to ^^C. Particularly, the Hoyle resonance (the second 0^ state) at 
7.65 MeV of ^^C plays an essential role in this triple-a process. At low temper- 
atures, say, T < 10^ K, however, the energy of three a particles cannot reach the 
Hoyle resonance. In this case, a direct fusion process of three a particles not through 
the Hoyle resonance becomes dominant. Thus, evaluation of the reaction rate of 
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is essential for understanding the formation of C. 
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Fig. 22. (left panel) Triple-a reaction rate as a function of temperature. The solid line represents 
the result of CDCC and the dotted line is the resuh of CDCC with the RSM. The dash-dotted 
line shows the reaction rate of NACREl — ^ (right panel) The solid and dashed lines show the 
reaction probability (aw)^ ^ of the nonresonant (ei = 38.2 keV) and resonant (ei = 92.0 keV) 
processes, respectively. The dotted line shows the result for the nonresonant process with the 
RSM in the CDCC calculation. 



In Ref. 1144]) . formulation of this ternary fusion process (TFP) based on CDCC 
was developed, and applied to the study of the nonresonant triple-a reaction. We 
describe resonant and nonresonant processes on the same footing. The left panel of 
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Fig. [22] shows the resulting reaction rate. The horizontal axis is temperature and 
the vertical axis is the order of the rate. The solid line is the new reaction rate 
calculated with CDCC, which is much larger than the rate of NACRE^^^' shown by 
the dash-dotted line. The difference is up to about 20 orders of magnitude around 
10^ K. We stress that it is shown in Ref. I144p that the approximate method to 
mimic the nonresonant triple-a process, the resonance shift method (RSM), proposed 
by NomotcP^ and used in NACRE, has no theoretical foundation. This can be 
understood clearly if one sees the reaction probability {av)ei^E, where ei is the a-a 
relative energy and E is the total energy of the three-a system. For clear discussion, 
we show the results of {(Tv)e^^E obtained with the coupled-channel effects switched 
off; the same model space as in the CDCC calculation converged is adopted. The 
dashed line in the right panel of Fig. [22] shows {av)e^^E for ei = 92.0 keV, i.e., the 
resonant capture probability. As for nonresonant capture process, we show by the 
solid line the result for ei = 38.2 keV calculated with CDCC, which has completely 
different energy dependence from that of the dashed line. On the other hand, in 
the RSM, the probability shown by the dotted line is used. One sees that it has a 
resonance peak at different energy from that of the Hoyle state [E = 387 keV). This 
is also the case with other values of ei. Therefore, one finds that the RSM implicitly 
assumes that there are infinite number of resonances around the Hoyle resonance, 
which is obviously inconsistent with the experimental information on ^^C. If we adopt 
in our CDCC calculation the assumption used in the RSM (and also in NACRE), 
we obtain the reaction rate shown by the dotted line in the left panel of Fig. [22] that 
agrees well with the rate of NACRE. 

It should be noted that the calculation of the triple-a reaction at low energies 
requires extremely large model space as described in Ref. I144p . We have confirmed 
a clear convergence of the reaction rate with respect to the model space of CDCC; 
the maximum value rmax of the a-a relative coordinate r is set to 5,000 fm. If r 
is truncated at a smaller value, say, 200 fm, we have a reaction rate much smaller 
than the solid line in the left panel of Fig. [22] This is also the case when a screened 
Coulomb interaction with a screening radius of a few tens of fm is used.l^^ It is also 
found that single-channel calculation never converges, and an adiabatic description 
of the three-particle system at low energies does not work at all. 

Our new reaction rate has been applied to several astrophysics s\^x^(XlQs^MEM^M^^^M^^M^ 
Some of them, Refs. I148P . 1149]) . llSOp . [TFT)) , show that our reaction rate is incompat- 
ible with observation, whereas others, Refs. ll52p . [T35]) . ll54p . claim that our rate is 
consistent with observation; a slight modification on our reaction rate is necessary 
in Ref. 1153] ). On the nuclear physics side, very recently, a new evaluation of the non- 
resonant triple-a reaction rate with the hyperspherical harmonic R-matrix method 
was performed by Nguyen and collab orators \^=^ They showed that the nonresonant 
contribution was indeed large, i.e., about 20 orders of magnitude larger than the rate 
of NACRE at 10^ K, which is qualitatively consistent with our finding. At temper- 
ature higher than 7 x 10^ K, however, their result agrees very well with the rate of 
NACRE. Further investigation on the nonresonant triple-a reaction rate, as well as 
experimental challenges to reveal the features of the low-energy three-Q continuum, 
will be very important for our understanding of the origin of ^^C. 
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§9. Applications of CDCC to nuclear engineering 

One of the advantages of CDCC is its fully quantum- mechanical aspect. There- 
fore, it is applicable to reactions at low incident energies, which are important for 
nuclear engineering. In this section, we show some results of applications of CDCC 
to studies in this field. 

9.1. Deuteron induced reactions on ^''^ Li 

In the international fusion material irradiation facility (IFMIF) project,^!^ deuteron 
induced reactions on ^'^Li are considered as one of the most promising reactions to 
generate high-intensity neutrons. Understanding of the ^''^Li(d, xn) processes, where 
X indicates all possible reaction products, at incident energies up to 50 MeV is 
highly important. Neutron spectra observed at forward angles show a broad peak 
with approximately half the incident energy.l^^ This suggests an importance of 
deuteron breakup processes, namel y, deute ron dissociation and proton stripping due 
to nuclear fields. In the past works ,'^^''-^^ these processes were treated using semi- 
classical models such as the modified intra-nuclear cascade (INC) modeP^ and the 
Serber model.'^^ Since the incident energy of interest here is relatively low, more 
sophisticated quantum mechanical approaches will be necessary for quantitative un- 
derstanding of the neutron production rate. 

The ^'^Li(d, xn) cross sections are expected to consist of the contributions from 
the deuteron elastic breakup and proton stripping processes. In Ref. I16ip we applied 
CDCC and the Glauber modeP^ to estimate the former and the latter, respectively, 
of the '^Li((i, xn) reaction at 40 MeV. As for the nucleon optical model potential for 
^Li, we took the parameter set determined in Ref. I162p through a systematic analysis 
of the nucleon elastic scattering on ^•'^Li. 



> 

S 



10" 



10" 



Neutron Production 
^ Proton Stripping (Glauber) 

10 Elastic Breakup (CDCC) 

cvap. + prceq. (Fitting) 




Hagiwara+2005 i — • — i '^-^.^ 
_i 1 I 1 1 1 1 I 1 1 1 i_l 



Neutron Production 
Proton Stripping (Glauber) 
Elastic Breakup (CDCC) 
evap. + preeq. (Fitting) 



10 20 30 

Neutron Energy (MeV) 



40 




40 60 80 

Neutron Angle (deg) 



Fig. 23. Comparison between the experimental data and the calculated result of d^ a / {dE\^dQ\^) at 
the neutron emission angle 10° (left panel) and the neutron angular distribution (right panel) 
of ^Li(d, xn) at 40 MeV. For the latter, the cross section is integrated over the neutron outgoing 
energies. The proton stripping and elastic breakup components are plotted by the dashed and 
dotted lines, respectively, and the sum of the evaporation and pre-equilibrium components is 
shown by the dash-dotted line. The solid line represents the total neutron production. 



In Fig. [23] we show the comparison between the calculated result and the ex- 
perimental data. The left panel shows the double differential breakup cross section 
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cf (J / {dEl^dQll) with respect to the energy E]^ and the angle i?^ of the outgoing neu- 
tron in the laboratory frame; the result corresponds to the neutron scattering angle 
at 10°. The angular distribution dcr/di?,^ is shown in the right panel. In each panel, 
the solid and dashed lines show, respectively, the contributions of the elastic breakup 
and the stripping. The dotted line is the sum of the evaporation and pre-equilibrium 

components, which are evaluated by an empirical method, i.e., the moving source 
model.nSJ 

The calculation reproduces the experimental data quite well. As expected, the 
deuteron breakup processes are strongly involved with the formation of the bump 
around 20 MeV of the outgoing neutron energy. The proton stripping process is more 
predominant than the elastic breakup process. It is worth noting that the strong 
contribution from the stripping process is also indicated by the DWBA calculation for 
{d,px) reactions on medium-heavy nuclei at 25.5 M.eV?^^ The relative contribution 
from the stripping process is reduced as the angle increases. It is found also that the 
deuteron breakup processes dominate neutron production at forward angles and are 
negligible at backward angles, while the evaporation and pre-equilibrium processes 
have a major contribution at backward angles. Consequently, the present analysis 
suggests that accurate description of the deuteron breakup reactions including both 
the elastic breakup and proton stripping processes is essential for reliable design of 
the neutron sources using the Li{d, xn) reaction. 

Although this model calculation reproduces successfully well the experimental 
data at small angles, we found some problems of the Glauber model calculation for 
the neutron stripping. The Glauber model calculation shows that the peak position 
in the emission spectra shifts to high energy as the emission angle increases, and fails 



to reproduce the experimental spectra at large angles; see Ref. I16ip . Furthermore, 
inclusive ( d, x p) spectra for heavy target nuclei are underestimated by the model 
calculatioipSiJ because the Glauber model cannot treat Coulomb breakup effects ac- 
curately. These may suggest a limitation of applying the Glauber model. Therefore, 
it is desirable to apply the eikonal reaction theory (ERT) described in ^ to the 
evaluation of the stripping process. This may eventually make it feasible to predict 
nucleon production from deuteron induced reactions. 

9.2. Neutron induced reactions on ^Li 

Lithium isotopes will be used as a tritium-breeding material in d-t fusion reac- 
tors, and accurate nuclear data are required for n- and p-induced reactions. Indeed, 
the international atomic energy agency (IAEA) is organizing a research coordination 
meeting to prepare nuclear data libraries for advanced fusion devices, FENDL-3,'^^^ 
and the maximum incident energy is set to 150 MeV to comply fully with the re- 
quirements for the IFMIF project. Lithium isotopes are some of the most impor- 
tant materials in these libraries. In fact, understanding of the breakup properties 
of lithium isotopes by neutron is crucial to determine the neutron energy spectra in 
blankets of fusion reactors. The tritium breeding ratio, nuclear heating distributions, 
and radiation damage of structural materials are affected by the n+Li reactions sig- 
nificantly. 

In spite of the importance of the n+^Li reaction mentioned above, experimental 
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data of the breakup processes leading to the ^Li continuum are extremely rare for 
the neutron energy region above 20 MeV. Furthermore, the statistical model, often 
used in evaluation of nuclear data for medium to heavy nuclei, cannot be applied 
to the ^Li(?7-,n') reactions. As mentioned in Refs. [T66]) and I167L the mechanisms 
leading to the three-body {n + d + a) and four-body {n + n + p + a) final states 
are particularly important. Therefore, more reliable theoretical calculations for the 
cross sections are highly desirable. 

In Ref. 1168] ). we performed a microscopic analysis of n+^Li scattering by means 
of CDCC. In this analysis, ^Li is described as a d + a system, and the interaction 
between n and ^Li is calculated by a folding model with the JLM effective inter ac- 
tion,l22' where the normalization of the imaginary part, A^, is optimized to reproduce 
the elastic cross sections. Details of the calculation are shown in Ref. 11681). 
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Fig. 24. Angular distribution of the elastic differential cross section of n+^'Li scattering for incident 
energies between 7.47 and 24.0 MeV. The solid and dashed lines correspond to the result with 
and without couplings to breakup states of ®Li, respectively. Experimental data are taken from 
Refs. n"5^ . ll70p . [17111 . The results are subsequently shifted downward by a factor of 10~^-10~^ 
from 8.96 MeV to 12.0 MeV on the left panel, and IQ-^-lO"" from 13.9 MeV to 24.0 MeV on 
the right panel, respectively. 

Figure [Ml shows the differential elastic cross sections of n+^Li for incident ener- 
gies between 7.47 and 24.0 MeV. One sees that the results of the CDCC calculation 
(the solid lines) are in good agreement with the experimental data. The dashed 
lines represent the results of a single-channel calculation, in which couplings to the 
breakup states are omitted. It is found that breakup effects shown by the differ- 
ence between the dashed and solid lines are significant to reproduce the angular 
distributions of the elastic scattering. For all incident energies, we take A^ = 0.1 to 
reproduce the data. It should be noted that the single-channel calculation cannot 
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reproduce the experimental data with any values of A^. The very small value of 
A^ obtained may indicate that most of the reaction channels that cause loss of the 
incident flux are treated explicitly in the CDCC calculation. 
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Fig. 25. The neutron spectra calculated by CDCC with the JLM interaction comparing with mea- 
sured data at selected angular points in the laboratory system. The experimental data are taken 
from Ref.nM]). 



In Fig. 1251 the calculated neutron spectra are compared with the experimental 
data at selected scattering angles in the laboratory frame and incident energies; 
A^ is fixed at 0.1. Components of the breakup states of ^S, ^D, ^D, and ^D are 
represented by the dash-dotted, dashed, dotted, and thin solid lines, respectively. 
These results are broadened by considering the finite resolution of the experimental 
apparatus.l^^ In each panel, three peaks of the experimental data correspond to 
the elastic scattering (right), the inelastic scattering to the 3"*" resonance (center), 
and the inelastic scattering to the 2+ resonance (left). The CDCC calculation gives 
a good agreement with experimental data in the high neutron energy region. On the 
other hand, in the low neutron energy region, which corresponds to highly exited 
states of ^Li, the calculated cross section considerably undershoot the experimental 
data. This indicates that experimental data contain contribution from the (n, 2n) 
process due to a four-body breakup reaction ^Li(n, nnp)a, as indicated in Ref. 1169] ). 
In the present calculation, the four-body breakup effects are treated only implicitly 
by using a complex optical potential. Estimate of these effects with four-body CDCC 
will be very important. 

Thus, CDCC with the JLM interaction is expected to be a powerful framework 
for the data evaluation of the ^Li(n, n') reactions. Once the JLM parameter is deter- 
mined by an analysis of elastic scattering, evaluation of the inelastic cross sections 
and neutron breakup spectra can be done with no free adjustable parameters. This 
is a very important feature of the present framework that enables a quantitative 
calculation of the cross sections for nuclear engineering studies. 
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§10. Summary 

This review has shown recent developments of the continuum discretized coupled- 
channels method (CDCC) and its applications to nuclear physics, cosmology and 
astrophysics, and nuclear engineering, after the previous review articles of Refs. [2l|) 
andESD- 

First, we have shown the theoretical foundation of CDCC. The primary ap- 
proximation in CDCC is the truncation of the angular momentum i between two 
constituents of the projectile. This ^-truncation changes the two-body potentials be- 
tween the target and the constituents to three-body potentials. The CDCC solution 
is the zeroth-order solution to the distorted Faddeev equations with the three-body 
potentials, the first-order correction to which is largely suppressed. The CDCC so- 
lution is thus a good solution to the three-body Schrodinger equation, when the 
maximum value of i is large enough. This theoretical statement has recently been 
confirmed numerically by directly comparing the CDCC and Faddeev solutions.'^ 
In CDCC calculations, it is necessary to confirm the convergence of the solution with 
respect to expanding the model space, so that the CDCC solution has no model space 
dependence and hence is a good solution to the three-body Schrodinger equation. 

As an underlying theory of CDCC, we have constructed a microscopic reaction 
theory for nucleus-nucleus scattering based on the multiple scattering theory. The 
input of the theory is either the t- or the ^r-matrix NN interaction instead of the 
realistic one. The former has much milder r dependence, so that the Glauber model 
becomes reliable for the systems in which the Coulomb interactions are weak, when 
either the t- or the gr-matrix is used as an input of the model. In the scattering of 
lighter projectiles from lighter targets at intermediate energies, effects induced by 
finite nucleus such as projectile breakup and target collective excitations are found 
to be small. Then the double-folding model becomes reliable for the scattering. 
Using the model, one can construct the microscopic optical potential with projec- 
tile and target densities calculated by fully microscopic structure theories such as 
antisymmetrized molecular dynamics (AMD), the Hartree-Fock (HF) method, and 
the Hartree-Fock-Bogoliubov (HFB) method. This fully microscopic framework has 
been applied to the scattering of stable nuclei and unstable neutron-rich Ne isotopes 
at intermediate energies with success in reproducing experimental data. The reliable 
microscopic optical potential can be used as an input of CDCC calculations, since the 
nonlocal potential can be localized with the Brieva-Rook (BR) method accurately. 
This microscopic version of CDCC is quite useful to analyze scattering of unstable 
nuclei systematically. 

Then we have extended CDCC in three directions, and developed the following 
new frameworks for breakup reactions: i) eikonal-CDCC (E-CDCC) for Coulomb- 
dominated breakup processes, ii) four-body CDCC for breakup of a three-body pro- 
jectile, and iii) the eikonal reaction theory (ERT) for inclusive breakup processes. 

E-CDCC treats both nuclear and Coulomb breakup very accurately and effi- 
ciently. The quantum-mechanical (QM) corrections, if necessary, can be made by 
replacing the partial scattering amplitudes corresponding to small impact parame- 
ters h with the solutions of a QM calculation, i.e., CDCC. E-CDCC calculations with 
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the QM corrections are much faster than CDCC calculations, with keeping the same 
accuracy as the latter. E-CDCC has thus become a standard method for describing 
breakup reactions in which the Coulomb breakup contribution is essential. Recently, 
inclusion of a relativistic Lienard-Wiechert potential as a Coulomb interaction was 
accomplished in E-CDCC, and dynamical relativistic effects generated by the poten- 
tial were found to give a sizable increase in the *B and ^^Be breakup cross sections 
at 250 MeV/nucleon. 

Four-body CDCC adopts pseudostates, which are obtained by diagonalizing the 
internal Hamiltonian of a three-body projectile, as discretized-continuum states of 
the projectile. The four-body wave function is expanded by the pseudostates together 
with the bound state(s); they are assumed to form a complete set in a space in 
which the reaction takes place. Additionally, a new method for obtaining a smooth 
breakup energy spectrum from discrete breakup cross sections given by CDCC was 
proposed. Four-body CDCC together with this new smoothing method based on the 
complex-scaling method has completed the description four-body breakup processes. 
Clear convergence with respect to expanding the model space is seen in both the 
elastic and the breakup cross section of the ^He scattering, and the results converged 
are consistent with experimental data. Four-body CDCC with the new smoothing 
method is indispensable for systematic studies of unstable nuclei that have a three- 
body structure. 

ERT is a framework that makes CDCC applicable to inclusive breakup processes 
such as neutron removal reactions, in which only the core nucleus of the projectile 
is detected. ERT gives separation of the scattering matrix into the contribution of 
each constituent of the projectile, without making the adiabatic approximation to the 
Coulomb potential. This is an essential point of ERT that eliminates the well-known 
shortcoming of the Glauber model, i.e., the divergence problem in Coulomb breakup. 
ERT has been successfully applied to the one-neutron removal from ^^Ne, and one- 
and two-neutron removal from ^He. With ERT, the accuracy of the Glauber model 
was systematically investigated for deuteron induced reactions at 200 MeV/nucleon. 
The Glauber model is good for light targets, but not for heavier targets, because of 
the Coulomb breakup contributions for the latter. 

CDCC, E-CDCC, four-body CDCC, ERT, and the microscopic version of CDCC 
have extensively and successfully been applied to studies of various reactions essen- 
tial in nuclear physics, cosmology and astrophysics, and nuclear engineering. These 
methods treat nuclear and Coulomb breakup processes nonperturbatively and nona- 
diabatically, and have no free adjustable parameters. These features are essential for 
quantitative studies such as the investigation of four-body breakup properties of ^He 
at the Coulomb barrier energy, determination of astrophysical S factors, evaluation 
of the contribution from the nonresonant triple-a process, the accurate estimation 
of neutron yields for the nuclear engineering design, and so forth. Validation of find- 
ings or conjectures obtained by nuclear structural studies through direct comparison 
with experimental cross sections is also highly important. The microscopic reaction 
theory plays definitely important roles in this subject, as we have confirmed ^^Ne to 
be a halo nucleus with large deformation. 

There are several interesting future works based on CDCC and the extended 
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theories. One is systematic applications of these theories to nuclear physics and the 
related fields. Another interesting work is to clarify the relation between CDCC and 
the dynamical eikonal approximation (DEA).'S2li'ElI Many nuclei can be described 
well by a core nucleus and extra nucleon(s). In the scattering of these nuclei, ex- 
citations of the core nucleus are important in general.'^^''^^ It is interesting to 
extend CDCC to treat the core excitations during the breakup. Use of other three- 
body structural models in four-body CDCC is also an important future work. The 
cluster-orbital shell model (COSM),'^^''^^ which has been applied to up to five- 
body systems, is one of the most promising models. Extension of four-body CDCC 
to five- and six-body CDCC with COSM will be a very attractive subject to explorer 
drip line nuclei. Finally, it will be very important to develop more accurate descrip- 
tion of transfer reactions at relatively low energies, by treating the rearrangement 
channels nonperturbatively. This will be essential to investigate excited states of 
stable and unstable nuclei, which are to be accessed via transfer reactions. These 
future works are important in further development of nuclear physics and the related 
fields. 
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